Small- x physics in perturbative QCD 



L.N. Lipatov ' 
Petersburg Nuclear Physics 
Institute, 

Gatchina, 188 350, St. Petersburg, Russia 



Abstract 

We review the parton model and the Regge approach to the QCD description of the 
deep-inelastic ep scattering at the small Bjorken variable x and demonstrate their rela- 
tion with the DGLAP and BFKL evolution equations. It is shown, that in the leading 
logarithmic approximation the gluon is reggeized and the pomeron is a compound state 
of two reggeized gluons. The conformal invariance of the BFKL pomeron in the impact 
parameter space is used to investigate the scattering amplitudes at high energies and 
fixed momentum transfers. The remarkable properties of the Schrodinger equation for 
compound states of an arbitrary number of reggeized gluons in the multi-colour QCD are 
reviewed. The gauge-invariant effective action describing the gluon-Reggeon interactions 
is constructed. The known next-to- leading corrections to the QCD pomeron are discussed. 
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1 Introduction 



Recent measurements of the structure functions for the deep-inelastic ep scattering at 
HERA discovered their dramatic rise at the region of small x ~ 10~ 4 [1]. In the framework 
of the Bjorken-Feynman parton model [2] this experimental result implies the correspond- 
ing growth of the parton distributions rii(x) inside the rapidly moving proton as functions 
of the decreasing parton momentum fraction x and the increasing photon virtuality Q 2 . 
In the framework of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation 
[3,4,5,6] the parton distributions grow at small x as a result of their (Revolution. In 
the framework of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation [7] this growth 
is a consequence of their x -evolution. Within the double-logarithmic accuracy these 
equations coincide and the increase of the structure functions at small x is related with 
the singularities of the anomalous dimensions for the corresponding twist- 2 operators at 
non-physical values j — > 1 of the Lorentz spin [2,8]. The existing experimental data on 
structure functions agree with the DGLAP dynamics provided that the evolution equa- 
tion in Q 2 is applied starting from rather small Q 2 = Ql [9]. The growth of the structure 
functions at small x can be also obtained with the use of the BFKL equation [10]. In this 
case a large uncertainty is related with the fact, that all next-to-leading corrections to 
this equation have not been calculated yet contrary to the case of the DGLAP equation 
where they are well known. 

The additional information on the dynamics of the deep-inelastic scattering at small 
x is extracted from the study of the final state particles. An especially clean footprint of 
the BFKL pomeron can be found in the processes with the inclusive production of jets 
[11]. The quark-gluon structure of the pomeron can be investigated at the hard diffractive 
scattering when the hadrons are produced in the virtual photon fragmentation region [12]. 
The deep-inelastic process with a large rapidity gap for the final particle momenta was 
discovered at HERA [13] . The various theoretical models for its interpretation were 
suggested [14]. The other high energy processes with the large rapidity gaps were widely 
discussed to discriminate the dynamics related with the soft and hard pomerons [15]. 

In this review we consider the theory of the BFKL pomeron. Because this theory 
is related closely with the parton description of the deep-inelastic scattering in QCD 
and with the Regge model, we remind below the basic ideas of these two traditional 
approaches to the high energy physics (a more comprehensive information can be found 
in Refs [16], [17]). In the next section the basic properties of the solution of the BFKL 
equation are discussed in the framework of the impact parameter representation . The 
gluodynamics is known to be a low energy limit of the super-string model of elementary 
particles which includes the quantum gravity. On the contrary one can expect that QCD 
at high energies could be described in terms of an effective field theory for string-like 
objects. In the third section it will be demonstrated, that in the Regge limit of large 
energies ^/s and fixed momentum transfers \f—i the gluon having the spin j — 1 at t — 
lies on the Regge trajectory j = j(t). Such reggeization property was assumed to be 
typical for hadrons which are extended objects. We consider here a simple effective field 
model in which the Feynman vertices coincide with the QCD reggeon-particle couplings. 
It is shown, that in the leading logarithmic approximation the pomeron is a compound 
state of two reggeized gluons. However, for restoring the S'-matrix unitarity one should 
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consider the contribution of the diagrams with an arbitrary number of the reggeized 
gluons in the t-channel. In the end of third section it will be demonstrated, that the 
equations for compound states of several reggeized gluons in the multi-colour QCD have 
remarkable properties: the conformal symmetry, the holomorphic factorization of their 
eigen functions and the existence of non-trivial integrals of motion in holomorphic and 
anti-holomorphic subspaces. The corresponding Hamiltonian turns out to be equivalent 
to the local Hamiltonian of the exactly solvable Heisenberg model with the spins being the 
generators of the conformal (Mobius) group. At high energies it is natural to reformulate 
QCD as an effective field theory for reggeized gluons. In the fourth section the gauge- 
invariant effective action for the interactions between the reggeized and usual gluons is 
constructed. The main results in the problem of finding next to leading corrections to the 
BFKL equation are reviewed in the fifth section. In Conclusion some unsolved problems 
are discussed. 



1.1 Parton model in QCD 

The deep-inelastic ep scattering at large electron momentum transfers q = p e — p e > and 
a fixed Bjorken variable x = ^ (Q 2 = —q 2 ) is a well investigated process for which the 
perturbative quantum chromodynamics (QCD) was traditionally and successfully applied. 
Before the QCD discovery the approximate scaling behaviour of the structure functions 
Wi t 2(x, Q 2 ) for the virtual photon-proton scattering at large Q 2 was derived in the frame- 
work of the Bjorken- Fey nman parton model [2]. In this model the transverse momenta 
kf- of partons inside the moving proton are assumed to be independent of Q 2 , the cross- 
section <jl for the longitudinally polarized virtual photon is zero and the cross-section 
ot for the transversally polarized photon is expressed in the impulse approximation as a 
sum of photon-quark cross-sections averaged with the distributions of quarks n q (x) and 
anti-quarks n q (x) in the proton [1]: 

a T = E^-^T^KW +ng(x)), (1) 

where e q is the quark charge measured in the units of the electron charge e and a = e 2 /47r. 
As a consequence of the charge and energy conservation the following sum rules for the 
parton distributions rii(x) inside the proton 

l = X) e <?/ dx(n q (x) -n q (x)) , 1 = ^2 I dxxrii(x) (2) 

are valid. 

In the renormalizable field theories the Bjorken scaling for structure functions is vi- 
olated due to the logarithmic terms (g 2 In Q 2 ) n appearing in the perturbation theory . 
The leading logarithmic approximation (LLA) for the structure functions corresponds to 
the sum of all such contributions [3]. The physical reason of the scaling violation is that 
because of divergencies one should introduce an ultraviolet cut-off A for the integrals over 
the transverse momenta k^ of partons being quanta of bare fields [4] : 

(?) 2 <A 2 . (3) 
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Nevertheless, the parton model representation (1) for the 7*p cross-section in terms of 
the quark distributions n q (x), riq(x) remains to be valid in LLA [4] if we identify the 
ultraviolet cut-off and the photon virtuality 

A 2 = Q 2 . (4) 

Moreover, one can express the inclusive probabilities rii(x) for finding the parton i with its 
momentum fraction x — k / \ inside the proton having the big momentum J? — > oo 
through the partonic wave functions ty n (ki,xi] k^x^ ■■■] k^,x n ) as follows 

^)=in^/nS^w 2 ^ 2 (E ^) so* - e *r) e s{x ri - Xi ) . (5) 



k 



Here the index enumerates n» partons of the type % in the state with n = J2k n k partons. 
In expression (5) and in the normalization condition for \l/ n : 

i=EnA/nS^ w * 2 (E *a - e (6) 

n i 'H- J r =l ) x r r r 

the most essential integral contributions at large A 2 = Q 2 in LLA correspond to the 
branching processes in which each virtual particle can decay into two others having signif- 
icantly bigger transverse momenta. The interference terms between the decay amplitudes 
describing the different branches are negligible in the physical light-cone gauge: 

n^Vf, = , nl = (7) 

for the gluon field v^(x) and therefore we can use the probabilistic picture to calculate 
various characteristics of the partonic cascade [4] . From the renormalizability of the theory 
the A-dependence of \^ n \ 2 is known: 

|tt„| 2 ~IlW > Eni = n » ( 8 ) 

i i 

where Zj are the renormalization constants for the corresponding parton wave functions. 
Using also the above probabilistic picture for essential contributions we obtain the follow- 
ing equation by differentiating the normalization condition (6) for the wave function \& n 
[4]: 



where 



Ni = f 1 dxnAx) (10) 
Jo 

is the averaged number of partons % in the proton. The quantity Pj d In Q 2 is the total 
probability of the parton decay during the "time" interval d In Q 2 and it depends on 
the effective charge g{Q 2 )- The normalization condition (6) should be valid for the wave 
function \& of each hadron, which due to relation (9) is compatible with the Callan- 
Symanzik equation [18] for the renormalization constants: 



4 



?—-z i = -P i Z i . (11) 



d In Q 

Analogously by differentiating the partonic expression (5) for rij(x) one can derive [4] 
the DGLAP equation [3-6]: 



d 



dlnQ' 



mix) = -P int (x) + J2 l l( ^rPk^)n k {x'), P t = Y. f dyP^ r (y) (12) 

7 Jx X X r JO 



governing the Q 2 -dependence of the parton distributions rii(x). In the right-hand side of 
this equation the first term describes the decrease of the number of partons of the type i 
as a result of their decay to other partons. The second term corresponds to the increase 
of rii(x) due to the decay of other partons into the states containing the partons of the 
type % . The decay probabilities Pk^i(xfx') and P k = J2i Jo dy Pk^i(y) are calculated in 
the form of the perturbative expansion over the running QCD coupling constant a(Q 2 ) 
at large Q 2 

A QCD 

Here A.q C d ~ 10 2 Mev is the fundamental QCD constant, N c = 3 is the rank of the gauge 
group SU(N C ) for QCD and rif is the number of the quarks with masses smaller than 
For example, in the lowest order corresponding to LLA one obtains [3-6]: 



xv, 
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(14) 



where (y) are proportional to the elementary inclusive probabilities for the transitions 
k — > i between the partons with definite helicities X k , A; (A = ±1 for gluons and A = ±1/2 
for quarks; A denotes the helicity A for the anti-quark). Other inclusive probabilities can 
be obtained with the use of the parity and charge conservation. In evolution equation (12) 
the infrared divergency at y = % — > 1 is cancelled. The dependence from interactions in 
the confinement region enters in rii(x) only through the initial condition at some Q 2 = Ql 
assumed to be sufficiently large: Ql ^> Aq CD . One can apply the evolution equation also 



5 



for distributions of quarks and gluons inside a virtual parton of the type k substituting 
7ij(x) — > D l k (x) and assuming, that 



D\{x) = 5 kl 5{x-l) 



(15) 



for Q 2 = Q\. The parton distributions in the proton can be expressed through D\ using 
the representation: 

mix) = eT^ fM Di£), (16) 

where the initial parton densities fk{(3) are not calculated in the framework of the per- 
turbation theory. The solution of the DGLAP equation for D l k (x) can be written in the 
following form: 



Dl(x) 



J a—ioo ZTTt \X / r 



dj fl\ 



(17) 



where in LLA 



Jc 



Ql q 2 4tt (3 2 [a(Q 2 ))> W) A MQVAJou) 



(18) 



The integration contour L = (o — too , a + ioo) in the complex plane of the t-channel 
angular momentum j is situated to the right of the point j '• = 1. The anomalous dimensions 
7 — lr{j) an d couplings Mj = M^(j) are determined from the secular equation: 



1 M i = 1 l{j)M k , 
where the matrix 7(7) is defined as follows: 

liU) = [ d y 

Jo 

and its matrix elements are given below [19]: 

1 



y 3 1 w l(y) - $kiY, w k(y) 



(19) 



(20) 



1/2^ _ N 2 C - 1 



7$(J) = 



- + 7 ~ x - 25, 
2 j(j + 1) 



> 7i/ 2 0') 



AT C 2 -1 



7i/ 2 0') 



1 , 7 1 1/2 o-) = ^, 7r 1/2 a) 



11 



7i 1 (j) = 2iV ( 



+ 



J + 2 
4j 2 + 4j - 2 



25'. 

6 ' (./- D.iU I DU I 2) 



iV c J'(J 2 -1)' 
2 

./'(./' • UU • 2) 
2 



^\j)=2N c 
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where 



(./ ')./(./ • Di./ • 2!' 
Sj = + 1) - ^(1) , 1>(x) = T'(x)/T(x) 



(21) 
(22) 
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and T(x) is the Euler T-function. Due to the charge and energy conservation the following 
sum rules are valid: 

7i/2(l) = , + 7l/2 (2) + 7 V2 (2) = , 

7l 1 (2)+ 7 r 1 (2) + 2n / ( 7l 1/2 (2) + 7 r 1/2 (2)) =0. (23) 

Note, that for the e + e~ annihilation the inclusive cross-section of the hadron produc- 
tion can be written in the framework of the parton model in terms of the fragmentation 
function D^(x) which is the inclusive probability to find the hadron h with the energy 
fraction x inside the parton %. These fragmentation functions are expressed through the 
distributions D\[x) of dressed partons k inside the bare partons i. In the leading loga- 
rithmic approximation the functions Df(x) and Df(x) coincide due to the Gribov-Lipatov 
relation [3]: 

= D*(x) . (24) 

Using the arguments based on the crossing symmetry and analyticity one can expect that 
these two functions are related also by the Drell-Levy-Yan relation 

Dftx) = ±xDl(-) , (25) 

x 

where the sign in its right hand side depend on the parton types % and k . In the gauge 
theory the point x — 1 is singular due to the Sudakov suppression of the quasi-elastic 
scattering. Nevertheless, this relation is fulfilled in a generalized form if in each order of 
the perturbation theory one would continue around x — 1 in D(x) the polynomials in 
ln(l — x) as the polynomials in ln(x — 1) [4, 20]. 

As it is seen from above formulas, for j — > 1 the anomalous dimensions 7 ^(j) = 
K1/2U) + 71/2 (j) and 7f U) = 7i 1 (j) + 7f 1 (j) have the pole singularities: 

7f (J) - 2 *Lzl , 7 * (J ) - 4 N c . (26) 

Therefore for x — > in the integrand of expression (17) for the distributions of gluons 
n g (x) and quarks n g (x) one obtains a saddle point at 



J 



\ 



X 



which leads after calculating the integral over j to the rapid growth of the total cross- 
section for the 7 *p-scattering [3,8]: 



°t ~ ^ exp yi6A^ c e In ^. (28) 

In next sections we shall show, that for sufficiently small x this asymptotics should be 
modified as follows 

11 4a 

a T ~ -= {-fB^KL ) UBFKL = _ Nc ln 2 (29) 

VQ x 71 
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in accordance with the BFKL equation [7]: 
d 



d In 



X 



Y n g (x, k±) = 2uj(k ± ) n g (x, k±) + / d k' ± K(k±, k' ± ) n g (x, k' ± ) (30) 



written for more general distributions n g (x,k±) depending on the longitudinal (x) and 
transverse (k±) gluon momenta. 

The next to leading corrections to the anomalous dimension matrix and to the splitting 
kernels of the DGLAP evolution equation in QCD were calculated in Refs [21] and [22]. It 
turns out [22] , that with taking into account these corrections the Gribov-Lipatov relation 
is violated. 

Another approach to the deep-inelastic ep scattering is based on the Wilson operator 
product expansion of two electromagnetic currents [19]. In this approach the 7*p scattering 
amplitude is presented at large Q 2 as a series over 1/x with the coefficients which are 
the moments of the structure functions. These moments can be expressed as a result 
of the operator product expansion through the product of the coefficient functions and 
the matrix elements of local operators O l between the hadron states. The coefficient 
functions depend on Q 2 and on an intermediate parameter /i being the normalization 
point for the running coupling constant a(fi 2 ). The matrix elements of O l also depend 
on the normalization point fi playing role of an ultraviolet cut-off in the corresponding 
Feynman integrals. In the product of the coefficient functions and the matrix elements 
the dependence from fi 2 is cancelled. In LLA we can chose fi 2 = Q 2 . Then one should 
calculate the coefficient functions only in the lowest order of perturbation theory. In this 
case for the polarized particles only the operators Q^ aiOL2 -- a 3 anc [ Q q a i a 2...a 3 j^y^g ^he 
Lorentz spin j and constructed from the quark fields ip{x) appear in the Wilson operator 
product expansion: 

7aia 2 ...oy = ^ x )^ D a2 D a3 ...D a ^( X ) , 7^ = 7^,7^ = 7*75 , (31) 

where Sym means the symmetrization over the indices cci, a 2 --.aj and the subtraction of 
traces. The covariant derivative D a equals 

D a = ^ + gv a (x) (32) 

where the anti-hermitian matrix v a is expressed in terms of the gluon fields belonging 
to an adjoint representation of the gauge group SU (N c ) 

v a (x)=t a v a a , [t a ,t b ]=f abc t c . (33) 

Here t a are anti-hermitian Gell-Mann matrices and f abc are the structure functions of the 
gauge group. The commutator of the covariant derivatives is proportional to the strength 
G^p of the gluon field: 

[D a , Dp] = g G a p , G a p = t a G a af3 = d a vp - dpv a + g [v a , vp] . (34) 

Operators (31) constructed from the quark fields have the canonical dimension d q = 
j + 2 and therefore their twist t — d — j equals t q = 2. There are other twist- 2 operators 
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Oq and Oq constructed from the gluon fields: 

(35) 

where tr implies the trace for the subsequent colour matrix. Under the renormalization 
the quark and gluon operators mix each with another and only certain linear combinations 
of them are multiplicatively renormalizable and diagonalize the evolution equations. It is 
convenient to use the light-cone gauge for the gluon field v. 

- n q + xp 2 

v. = v a n a = 0, n = — , n a = 0. (36) 

2pq 

In this gauge the light-cone components of the operators O aia2 ' aj appearing in the opera- 
tor product expansion are the simplest examples of the quasi-partonic operators [23] which 
have two properties: a) they are universal for all possible interactions, which means in par- 
ticular that the quasi-partonic operators do not depend on the QCD coupling constant g ; 
b) their matrix element between hadronic states are expressed through the integrals from 
their partonic matrix elements on the mass shell averaged with the partonic correlators. 

Let us consider the twist- 2 operators for a general case, when their matrix element is 
calculated between the initial and final hadron states with different longitudinal momenta 
[24]. Then we can construct non-equivalent operators as bilinear combinations of fields 
with a different number of derivatives: 

0\ .... (x) = n ai n a *...n?0 I q aia *" a * = £ c k ^(id .) V ,71 = 7^,7^ = 7«75 , 

oh....= E d * <^{id )\idy- k nf v a ^ , = sf , = ef (37) 
k=i 

where d means, that the derivative acts at the function situated in the left hand side 
from the differential operator. The coefficients c r , d r can be fixed up to a common factor: 

(_i)fe 

° k = fc!(jfe + l)!(j- fc-l)!(j-fc)! ' dk = (k - l)\(k + - k + - k - 1)\ ' (38) 

from the requirement, that the operators O q and Oq are irreducible under the conformal 
transformations [25]. It is enough to verify their covariant properties under the inversion 
x a — > x' a = x a /x 2 with the simultaneous substitution of the fields: 

ip(x) — > ip'(x') = x 2 xip(x) , v^x) — > v'Jjc') = x 2 A llu (x) v v {x) , (39) 

where 

A^(x) = ^-2^. (40) 

Ob 

The free Green functions for the fields ip(x) and F a p(x) = d a vp{x) — dpv a {x): 

x-y 



(0 T^{x)^{y) 0) 



(x - y y ' 
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(0\TF, u (x)F a p(y)\0) 



A^aO - y)Kp{x -y)- A^(x - y)A ua {x - y) 



(41) 



(x - y y 

are transformed properly under the inversion. Analogously the matrix elements of the 
twist-2 operators (37) with coefficients (38) can be found in the free theory: 



(0 



(0 



TF^(x)0 I a _(z)F^(y) 



0) 



0) 



x 



y. - z. 



(x-z) 2 {y-zf 



x — z 



'7. 



X 



y. - z. 



(x-z) 2 (y-zf 



{x-zY '■ (z-y)*' 

3 ~ 2 ^ (x. ~ zf (y. - zf 
(x — z) 6 (y — zf 



(42) 

and are transformed also correctly [25]. 

Conformally-covariant operators have simple multiplicative properties under the renor- 
malization in LLA in accordance with the fact, that their matrix elements between the 
parton states in the momentum representation are proportional to the Gegenbauer poly- 
nomials which diagonalize partly the Brodsky-Lepage evolution equations generalizing the 
DGL AP equations for a non-zero momentum transfer [24] . These equations are invariant 
under a sub-group of the conformal group which includes the transformation 



x a - fit , < = (43) 

1 + 2 e x a n a 

where n a is the above light-cone vector and e is the group parameter. One can verify, 
that the gauge condition v a n a = is compatible with this conformal subgroup. 

For the simplest supersymmetric Yang-Mills model containing apart from the gluon 
also a gluino being the Majorana particle belonging to the adjoint representation of the 
gauge group the anomalous dimensions can be easily obtained from the above QCD ex- 
pressions. If we introduce the quantities 

7^0')=7AfO')±7A*0') (44) 
related to the anomalous dimensions of the operators O v and A the corresponding 
matrices can be parametrized by two functions 7i(j) and 72 (j): 



7f0') 



7f 0') 
7f0') 



7i(i) 
2j + l 




7|-(j) 



U) i 9 g ~(j) 



72 (J) 
2j + l 



J -I J -I 
J +2 j + 2 



7iC? + l) 
2j + l 



J + 2 -j + 1 
-J-2 3-1 



(45) 



This parametrization is a consequence of the fact [20], that there are four irreducible 
representations of the super-conformal group for the corresponding twist-2 operators with 
the degenerate anomalous dimensions proportional to 71 and 72 : 



7i U) = N c 



33 



, 7 2 (j) = N c 



2 Sj+i ~ o - 



(46) 
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In particular from the above formulas we obtain the Dokshitzer relation [6] for the 
anomalous dimension matrix elements 

l q q ti) + l 9 q ti) = l q 9 ti) + l 9 9 ti) (47) 

corresponding to a particular case of the Shmushkevich rule stating, that the probability 
of a process induced by a particle belonging to a group representation is equal for any 
member of this representation provided that all final states related by the symmetry 
transformation are taken into account. Because 7|(j), 7f(j) and ~f*{j) coincide in QED 
[4] and in QCD [5, 6] up to simple colour factors, from the Dokshitzer relation (47) we 
conclude, that the pure Yang-Mills anomalous dimension 7|(j) (and the corresponding 
splitting kernel w 9 g (x)) could be obtained from the QED results. 

In the next-to-leading approximation the conformal invariance is violated due to the 
conformal anomaly related with the non- vanishing (5 -function [26]. Nevertheless the 
solution of the Brodsky-Lepage equations can be found if one would take into account the 
Ward identities of the conformal group [27]. 

The quasi-partonic operators [23] of an arbitrary twist seem to be a natural general- 
ization of the above twist- 2 operators related with the probabilistic picture. Their matrix 
elements can be presented only in terms of the parton correlators being the integrals from 
the product of initial and final state wave functions with a different number of partons. 
They do not mix with other operators under the renormalization. The evolution equa- 
tions for the quasi-partonic operators have the form of the Schrodinger equations with 
a pair-wise interaction. The pair Hamiltonians are proportional to the splitting kernels 
for twist-2 operators with generally non-zero quantum numbers including colours, fla- 
vors, spins and momentum transfers. All these kernels are calculated in LLA [23]. They 
have a number of remarkable properties including their conformal invariance. The evolu- 
tion equations are simplified at the small-x region, which gives a possibility to find their 
solutions in the double- logarithmic approximation [28]. 

The higher twist contributions are important for finding screening corrections respon- 
sible for restoring the S'-matrix unitarity. This problem will be discussed in this review 
later in more details. Here we want to remark only, that one can calculate correlators 
of parton densities within a certain accuracy using the probabilistic picture [29]. Let us 
introduce the exclusive multi-parton densities using the following definition: 



/(4 1) ,4 1) .-.^ ) ;4 2) ...< ) ) — 7T~~\ /il 4)7tv, l^l 2 (2-WE^V(i - E-S r) ) • 

(48) 



Ml 

Then the generating functional for these quantities 

m = e / nn^ (r) <pM) (49) 

J r i 

satisfies an evolution equation in Q 2 which can be obtained from the above probabilistic 
picture. If we consider only the gluodynamics and take into account the most singular 
part of the corresponding splitting kernel in (14) by putting 

4 N r f 1 AN r dx 



gi \ - " 'c 



tf x (1 — x) Jo i(l-i) 



f 1 4 7V C dx 

L i7T=i) • (50) 
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these equations can be easily solved: 



T( ) = r°° dl exp (-/) w exp (-Q (1 - exp(-wQ) 1 
W) J-ioo 2ni 1 - ^ (1 - exp(-<)) J °° f exp(-xZ) ' 1 ] 

The inclusive multi-gluon densities are obtained from this expression by the change 
of its argument: ip — > 1 + x an d expanding the obtained generating functional in the 
series over x( x )- I 11 particular for x — we obtain: 1(1) = 1, which corresponds to the 
normalization condition for the partonic \l/-function. The generating functional 1(1 + x) 
for the inclusive partonic correlators can be used for finding the screening corrections at 
the quasi-Regge region x — > 0. A more general method taking into account the possibility 
of the parton merging was developed in ref . [30] . 

Below we survey the basic results of the Regge approach used in the BFKL theory. 

1.2 Regge theory 

The elastic scattering amplitude A(s,t,u) for the process a + b — > a' + b' is an analytic 
function of three invariants [17]: 

S = (Pa + Pbf, U =(Pa~ Pb') 2 , t= (Pa- Pa'f \ ( 52 ) 

related as follows: 

s + u + t = 4m 2 (53) 

for the case of spinless particles with an equal mass m. 

The function A(s,t,u) describe simultaneously three channels. In the s, u and t 
-channels one has 

s > 4m 2 , u < 0, t < 0; 

u > Am 2 , s < 0, t < (54) 

and 

t > 4m 2 , u < 0, s < (55) 

correspondingly. For example, the scattering amplitude in the s-channel (where y/s is 
the c.m.energy of colliding particles) is obtained as a boundary value of A(s,t,u) on the 
upper side of the cut at s > 4m 2 in the complex s-plane. 

The Regge kinematics of the scattered particles in the s-channel corresponds to the 
following region: 

s — — u » m 2 a* —t = if 2 , (56) 

where if is the momentum transfer in the cm. system (q — p a — p a ,). From the S'-matrix 
unitarity one can obtain the optical theorem for the total cross-section: 

a tot (s) = -Im s A(s,0), (57) 

s 

where Im s A(s, t) is the imaginary part of the scattering amplitude in the s-channel. The 
whole amplitude is expressed through its imaginary parts in the s and u channels with 
the use of the dispersion relation : 

1 f°° 1 1 f°° 1 
A(s, t) = - ds' ImjA + - / du Ira,, A, (58) 

71 J4m2 S' - S-16 71 Jim? V,' - U - 16 
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where we omitted the possible subtraction terms. 

In the Regge kinematics (56) the essential s-channel angular momenta I = pp (p is 
the impact parameter and p is the cm. momentum) are large. Therefore one can write 
the angular momentum expansion of A in the form of its Fourier transformation : 

A(s, t) = -2is J d 2 p [S{s, p) - 1] ef$~$ . (59) 

The s-channel partial wave S(s, p) can be considered in some models as a two-dimensional 
S- matrix which is parametrized by the impact parameter ~j? (see for example [31] ). 
Usually it is written in terms of the eikonal phase 5(s,p): 

S(s,p)=e i5M . (60) 

In accordance with the s-channel unitarity one obtains Im 5 > 0. Because essential values 
of p can not grow more rapidly than p max = cln(s) the Froissart theorem is valid: 

a tot <A7rc 2 ln 2 (s). (61) 

If the scattering amplitude is pure imaginary at high energies one can derive from disper- 
sion relations (58) the Pomeranchuck theorem for the particle-particle and particle-anti- 
particle total cross sections: 

of = of. (62) 

In a more general case when its real part is as big as possible this theorem is modified as 
follows 

of /of = 1, Re(A) oc Im{A) oc sln 2 {s). (63) 

In the Regge model the asymptotics of the elastic scattering amplitude in region (56) 
has the following factorized form 

Ms,t) = ^Mtf Mt) 9m9m- (64) 

V 

Here gi^if) are the Reggeon couplings with external particles, is the signature factor 
(for the signature p = ±1): 

. COS 7TJ 4- p 

q = i : — (65) 

J Sin 7TJ 

and j p (t) is the Regge trajectory assumed to be linear: 

j p (t) = f Q + a' p t, (66) 

where jg and a' are the reggeon intercept and slope correspondingly. 

For the case, when the trajectory j p (t) passes through the physical value j = n (differ- 
ent for two signatures p = ±1) corresponding to the integer (or half integer) spin a = n 
for an intermediate state, the amplitude takes the form: 

s n 

A(s,t)<x-—-, (67) 
z z n 
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where t n is the squared mass of the compound state lying on the Regge trajectory [17]. Ex- 
perimentally all hadrons constructed from light quarks belong to the Regge families with 
almost linear trajectories and an universal slope a' ~ lGev ~ 2 . As it will be demonstrated 
below, in the perturbative QCD the gluons and quarks are reggeized, which means, that 
they lie on the corresponding Regge trajectories [7]. 

The Regge asymptotics seems to be natural [17] from the decomposition of the scatter- 
ing amplitude in the sum of contributions of various angular momenta j in the t-channel 
where the scattering angle 9 t is related with s as follows 

Z = COS l + 2s ( 68 ) 
t-Am 2 

For its symmetric and anti-symmetric parts 

A(s, t) = A + (s, t) + A~(s, t), A ± (-s, t) = ± A ± (s, t), (69) 

this decomposition continued analytically to big s and fixed negative t takes the simple 
form: 

A p (s,t)= ^-e 3 s j m (70) 

Jo— zoo ZTli 

and satisfies asymptotically the dispersion relation (58). The functions <fij(t) are propor- 
tional to the t-channel partial waves fj(t): 

Mt) = c 3 p{Am 2 - t)~' fj(t), c 3 = 16nH> f^jfrj ( 71 ) 

and are real in the physical region of the s-channel. 

Thus, for the imaginary parts of the signatured amplitudes A p one obtains the simple 
formulae corresponding to the Mellin transformations: 

Im.A p (8,t)= r* 00 (72) 

The inverse Mellin transformations 

0P = ^ d£e- j t Im s A p (s,t), £ = ln(s) (73) 

are simplified versions of the Gribov-Froissart formulae for ff(t) [17]. 

The t-channel elastic unitarity condition for partial waves analytically continued to 
complex j for t > 4m 2 takes the form 

+ = cJ i {t _ Am yH t -h m + ie) m _ ie) . ( 7 4) 

Using the analytic continuation of this relation to complex t one can express the 
amplitude ^ on the physical sheet of the t-plane through its value in the same point on 
the second sheet and obtain for the t-channel partial wave fj(t) the Regge poles [17]. 
But the t-channel partial wave could also have fixed square-root singularities of the type 
<fij(t) = a(t) + b(t)y/j — jo as it takes place in renormalizable field theories including QCD. 
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Total cross-sections for hadron-hadron interactions are approximately constant at high 
energies (up to possible logarithmic terms). To reproduce such behaviour in the Regge 
model a special reggeon is introduced. It is called the Pomeranchuck pole or the Pomeron. 
The Pomeron is compatible with the Pomeranchuck theorem (62) because it has the 
positive signature and the vacuum quantum numbers. Its trajectory is assumed to be 
close to 1: 

j(t) = l + u(t) , uj = A + a't, (75) 

where A^ 0.08 and a' ~ 0.3 Gev~ 2 [32]. It means, that the signature factor approximately 
equals i and the real part of A(s, i) is small in the agreement with experimental data. 
The above quantities A and a' are the bare parameters of the so called soft Pomeron. The 
BFKL Pomeron has a big value of A. 

In the j-plane there should be other moving singularities of 4>^{t) - the Mandelstam 
cuts arising as a result of multi-Pomeron contributions to the t-channel unitarity equa- 
tions. V.Gribov constructed the reggeon diagram technique in which all possible Pomeron 
interactions are taken into account [33] . In the reggeon field theory the parameters of the 
initial Lagrangian are renormalized. The simple model for taking into account the con- 
tributions from the multi-Pomeron exchanges is based on the assumption that the phase 
in the eikonal representation can be calculated by the Fourier transformation of the am- 
plitude written in the Regge form : 

S(s, P) = yj0ji * 9 2 {t) (76) 

In this case the resulting amplitude satisfies the Froissart requirements. The analogous 
unitarization procedure can be used also in other cases when the scattering amplitude 
obtained in some approximation grows more rapidly than any power of In (s) (see [34]). 

The high energy theorems do not forbid the existence of another Regge pole with 
the vacuum flavour quantum numbers - the Odderon which has the negative signature 
and the negative charge parity [35]. It could be situated also near j = 1, which would 
lead to a large real part of scattering amplitudes at high energies and to a significant 
difference between proton-proton and proton-antiproton interactions. Such singularity 
appears in the perturbative QCD simultaneously with the Pomeranchuck singularity [36] 
and therefore the discovery of the Odderon effects would be very important. 

Due to the optical theorem (57) the total cross-section is proportional to the imaginary 
part of the scattering amplitude, for which the Regge asymptotics is usually postulated. 
Therefore it is natural to ask what production processes are most probable at high ener- 
gies. In the dual resonance model non-vacuum reggeons in the t-channel are obtained as 
a result of summing over the s channel contributions from resonances with growing spins 
and masses and the Pomeron is dual to the s-channel background. In multi-peripheral 
models all reggeons are constructed as t-channel compound states of partons and from 
the s-channel point of view they describe the processes of the multi-particle production. 
The produced particles have the multi-peripheral kinematics: their transverse momenta 
kf- are fixed and their longitudinal momenta are ordered in such way, that the invariant 
squared energies of neighbouring particles s« = (ki + k^i) 2 are also fixed. Therefore the 
average number of particles in the multi-peripheral models grows as ln{s). 
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In QCD the situation looks to be simpler than in the multi-peripheral models. Namely, 
here the average transverse momenta of produced particles grow slowly with energy and 
therefore due to the asymptotic freedom the effective coupling constant decreases. As a 
result in QCD (at least in the perturbation theory) the most essential contribution to the 
total cross-section in the reaction 

a + b — > a' + di + d 2 + ... + d n + b', a' = d , b' = d n+1 (77) 

arises from the multi-Regge kinematics for produced particle momenta: 

s > Si = (ki + /Ci_i) 2 > m 2 , -U =q! 2 ~ m 2 , (aJ) 2 ~ m 2 . (78) 

There is the following constraint for energy invariants Sj resulting from the reality condi- 
tion k 2 = m 2 for the final state particles 

n+l n ¥ 

n^ = *ii[™ 2 +(^) 2 ]- (79) 
i=i i=i 

The momentum transfers qi in the crossing channels t« are expressed through external 
particle momenta as follows: 



Qi = Pa- Pa> ~J2 k r- ( 80 ) 
r=l 

In terms of the Sudakov parameters for produced particles 

k = (3iP a + ctiPb + k^~, (kl,p ajb ) = 0, k 2 = satifii - (kf) 2 = m 2 (81) 
the multi-Regge kinematics looks especially simple: 



2 2 

m m 



1 > /3i > fa > ••• > Pn > — , — < Oil < a 2 < ... < On < 1, (82) 

s s 

saS ~ (A?) 2 ~ m 2 . (83) 

The decomposition of the momentum transfers g« in terms of the Sudakov parameters for 
produced particles is also simplified in the multi-Regge region: 

Qi = PiPa ~ ai-iPb + qt, Qi - -{qt? ■ (84) 

Thus, momentum transfers with a good accuracy are transverse vectors. Finally the 
energy invariants also can be expressed through the Sudakov variables: 

Si = sfc-ioti (85) 

and therefore the produced particle momenta in the multi-Regge kinematics are essentially 
longitudinal. 

Now we review shortly the theoretical description of the multi-Regge processes in the 
framework of the Regge model. It is natural to generalize the formulae for the elastic 
processes in the form: 

iW 2 = s^s^\..si { tr ] l(qu 92, -<M + i)- (86) 
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Here the function 7(^1, ...q n+ i) is not real and should contain something similar to signa- 
ture factors for the elastic amplitude to satisfy analytic properties in the direct channels 
Si. These properties are significantly simplified in the multi-Regge regime, as it was 
shown by K. Ter-Martirosyan, H. Stapp, A. White, J. Bartels and others. Namely, the 
inelastic amplitude has only physical singularities in the channels where the real produc- 
tion of intermediate particles is possible. There can be simultaneous singularities only in 
non-overlapping channels. 

Let us consider the amplitude A 2 _>3 for the single particle production in the multi- 
Regge kinematics. The signatures in channels t\ and t 2 are assumed to be equal to p\ and 
P2 correspondingly. The signatured amplitude has two contributions satisfying the double 
dispersion representation in the non-overlapping channels S\,s and s 2 , s correspondingly 
and can be written with the use of the double Watson-Sommerfeld transformation as 
follows: 



3 ~J 2ni2ni\- ^' 1 ^ 2 ^ 2 ^ 2 2 { Hi s ^hh\ > v*') 



1 2 m 

where two partial waves 4>j[j 2 are real functions. They are expressed through the spectral 
functions p(si,s) = Im Sl Im s A2^z and p(s 2 ,s) = Im S2 Im s A2^z with the inverse Mellin 
transformations: 

^nn = / dZwe-'^i* / dfe"^ p(a li2 , s), 6,2 = Zn(a li2 ), £ = ln(s). 

JO "^£l,2 

Using the unitarity conditions in the direct channels the spectral functions can be pre- 
sented in a non-linear way again through inelastic amplitudes, which leads to a set of 
non-linear equations. Such approach is convenient for finding the high energy asymp- 
totics in non-abelian gauge theories [7]. 

1 2 

In the framework of the Regge theory (j)j[j 2 are given in the factorized form: 

<l&=9(ti)- ^-TprTWfahX)-. 1 -Tn^)- (89) 

Ji-J Pl {ti) 32-3vi\ t V 

Here V is the reggeon-reggeon-particle vertex depending on the usual invariants t 12 and 
the produced particle transverse momentum k± expressed through s 1 , s 2 with the use of 
the reality condition 

kl 2 = ^ - m 2 . (90) 

s 

Thus, the production amplitude in the Regge model have the form 

A^ = fr< tl) fr™i(q 1 ,< b ), (91) 

where the real and imaginary parts of 7 are expressed in terms of the vertices g and T. 

In a general case of the amplitude (86) for n-particle production in the multi-Regge 
kinematics one should introduce many partial waves 4 >t j l ...j n+l describing different disper- 
sion contributions from all non-overlapping channels. 
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2 BFKL pomeron in the impact parameter space 



In this section we remind comparatively old results concerning the description of the 
BFKL Pomeron in the impact parameter representation [37]. 

The asymptotic behaviour of scattering amplitudes in the Born approximation is gov- 
erned by the spin a of the particle exchanged in the crossing channel: 

A Born ~ s CT (92) 

and the Regge asymptotics is a generalization of this rule to continuous values of the 
spin: a — > j = j(t). In higher orders of the perturbation theory the scattering amplitude 
behaves as s n (apart from possible logarithmic terms), where the power n — 1 +J2i (°i — 1) 
grows linearly with the spins <7j of the particles in the t -channel intermediate states. 

In QCD the gluon spin a is 1 and therefore here the most important high energy 
processes are caused by the gluon exchanges. For example, the Born amplitude for the 
parton-parton scattering is [7] 

A(s, t) = 2s g5 Xa;Xa , T% A - t g5 XbtXb , T%, B , (93) 

where \ are helicities of the initial and final particles; A, A', B, B' are their colour indices 
and T?j are colour group generators in the corresponding representation. The s-channel 
helicity for each colliding particle is conserved because the virtual gluon in the t-channel 
for small q interacts with the total colour charge Q c commuting with space-time trans- 
formations. 



2.1 Impact factors 

Let us consider now the high energy amplitude for the colorless particle scattering de- 
scribed by the Feynman diagrams containing only two intermediate gluons with momenta 
k and q — k in the t-channel. With a good accuracy we can neglect the longitudinal 
momenta in their propagators (cf. (84)): 

k 2 ~ k\ , (q - kf ~ (q - k) 2 ± . (94) 

The polarization matrix for each gluon can be simplified at large energies s = {p a +Pb) 2 ^> 
m 2 as follows 

5^ = 5{? v + 5T ~ Sr = p * p " b+p * p b . (95) 

" ~ " PaPb 

The projector to the longitudinal subspace 5^ u can combine the large initial momenta 
p a ,Pb in a big scalar product s/2 = p a Pb- Moreover, if the indices ji and v belong to the 
blobs with incoming particles a and b correspondingly, then with a good accuracy we have 

5 „u _^ PbK^ / 96 x 

PaPb 

From the point of view of the t-channel unitarity for the partial waves fj(t) with 
complex j this substitution has a rather simple interpretation. The most important 
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contribution in the t-channel appears from the nonsense intermediate state leading to a 
pole singularity of fj(t) at j = 1. For this nonphysical state the projection of the total spin 
S^to the relative momentum of gluons in the cm. system equals 2 which corresponds to 
the opposite sign of their helicities. It means, that here the complex polarization vectors 
ei and €2 of the gluons coincide each with another: e\ = €2 = e. The projector to the 
nonsense state is e^e^ e*e*/ (note, that e 2 = 0). Due to the Lorentz and gauge conditions 
the vectors e are orthogonal to both gluon momenta k and q — k. Because these momenta 
are almost transverse at high energies (see (94)), after the analytic continuation from 
the t-channel to the s-channel the above projector should be proportional to the product 
PbPb PaPa i n accordance with (96). 

Using eq.(96) and introducing the Sudakov parameters 

a — — —s a / s , (3 — — — = Sb/s , k = k _|_ , d 4 k = d 2 k a< ^ b (97) 

PaPb PaPb 2 \s\ 

for the virtual gluon momenta k ,q — k one obtains for the asymptotic contribution of the 
diagrams with two gluon exchanges the following factorized expression: 

A(s, t) = 2t \s\ I [d 2 k^ 1 $ a (t, If-t) $ 6 (t, t-t), (98) 

k {q-k) 2 

corresponding to the impact-factor representation [34]. Here the sum over the colour 
indices is implied and the factor 1/2! in front of the integral is related with the gluon 
identity and compensates the double number of Feynman diagrams appeared due to our 
definition of the impact factors $ a,b as integrals over the energy invariants s a ^b from the 
photon-particle amplitudes fffl: 

t-t) = r p^ p -k?k ff( 8aJbi t, t-t) . (99) 

J-00 {ZTC) z t s s 

The impact factors describe the inner structure of colliding particles. For large k they 
are proportional to the number of partons iVj weighted with their colour group Casimir 
operators. 

Note, that we neglected the longitudinal momenta in gluon propagators in representa- 
tion (98) in accordance with (94) because in the essential integration region for the impact 
factors we have 

s a>b ~ m 2 , ~t 2 ~ (lf-~fc*) 2 ~ m 2 (100) 

and therefore k 2 = k\. The ultraviolet convergency of the integrals over s a , b 

follows from the fact, that due to the Ward identities k^f^ = (q — k) v f^ v = for the 
scattering amplitudes / = f a ' b one can perform the substitution: 



jCpSU^ klti_JYL f;t (101) 



ja,b 

S S Sn b Sn b 



Indeed, it has been assumed above that the amplitudes / do not contain pure gluonic 
intermediate states in the t-channel and therefore f a,b s a>b at high energies, which 
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leads to the rapid convergency of the integrals over s a ^. Therefore one can enclose the 
integration contours around the right cuts of f a ' b in the complex s a) b planes: 

if-t) = r Imsa f {s t, (102) 

Jth 2tt z s a , b s a , b 

after calculating its discontinuity. From this representation of <£> we conclude that the 
impact factors are real functions of k,lf—k, vanishing for small | k | and \lf — k \ 
in the case of the colorless particle scattering (e.g. photon-photon collisions), which is 
a consequence of the gauge invariance and of the absence of infrared divergencies in the 
integral over s ab for small | k | and |~if — k \. 

From the physical point of view the infrared stability follows from the fact, that for 
small | k | the virtual gluon interacts with the total colour charge of scattered particles. If 
this charge is zero, then the dipole and (generally) multipole interactions are proportional 
to powers of k . In particular the total cross-section for the photon-photon interactions 
does not contain any infrared divergency in the integral over k . The various impact- 
factors for the real and virtual photons are calculated [7, 34], which in particular gives 
us a possibility to estimate the hadron-hadron cross-sections using ideas of the QCD 
sum rules [38]. For hadrons the impact factors can be extracted phenomenologically for 
example from the inclusive lepton-hadron scattering because apart from the colour and 
charge factors these quantities are the same for the virtual gluons and photons. 

2.2 Mobius invariance of the BFKL pomeron 

It is convenient to present eq.(98) in the form of the Mellin transformation (70) 

A(8,t)=i\s\ J^^f u (^),t = -^ ( 103 ) 

and to pass to the impact parameter representation (cf. (59)) performing the Fourier 
transformation [38]: 

J r=l,2 y^ n ) 

(104) 

The quantity f u (pi, p|; p^,p^) can be considered as a four-point Green function: 

Un, ti; tf, p?) = (o |0(pi) 0(p 2 ) <KPrO 0(p 2 ')I °) ( 105 ) 

where the field <p(p) describes the (reggeized) gluons. In the Born approximation f w is 
proportional to the product of free Green functions: 

_ 4tt 2 

/"(Pi, P2', P?, P?) = ln \Pll' I ln \P 2 2' I » Pit = Pi ~Pk, (106) 

UJ 

related with the gluon propagators 

4i = -/^exp(<t^)ln|p|. (107) 



20 



The functions <3> a ' 6 (pi, pi, if) are related with the impact factors by the Fourier transfor- 
mations: 

<$> a '\pl, P2, ~f) = J d 2 k <!> a ' b (t, -f-lc) exp(^pt) exp(i(^ - t)p%) . (108) 

The vanishing of $ a,fe ( k ,~f — k) (102) at k — > or at fc — > if is equivalent to the 
following sum rules for $ a ' b (pi, p|, if): 

J $ a ' b (pt, pi, If) d 2 Pl = J $ a ' b (pt, ti, t) d 2 P2 = . (109) 

Therefore the expression for f w (q 2 ) is not changed if we add to fui{p~\,p~2',p~f,p~^) an 
arbitrary function which does not depend on pi, pi, pf or p^t. We can use this freedom 
related with the gauge invariance to modify /°(pi, pi; pjsP?) in the following way: 



2vr s 



In 



Pll'P22' 



In 



PnPyi 



, Pik = Pi — Pk 



(110) 



Pl2'Pl'2 

This expression is unique in comparison with all other physically equivalent expressions 
for /° because it depends only on two independent anharmonic rations of the vectors pi, 
p|, pT' and p^ which can be chosen as follows 



a 



Pll'P22' 



Pl2'Pl'2 



Pll'P22' 



PyiPv-2' 



, 7 







Pl2'Pl'2 



Pl2Pl'2' 



Therefore /° is invariant under the conformal (Mobius) transformations: 

a Pk + & Pfc 

Pfc — y 

cpk + dp k 

for arbitrary complex a, 6, c and d provided that we use the complex coordinates 



(111) 



(112) 



(113) 



Pk = x k + iy k , p k = x k -iy k 

for all two-dimensional vectors p~k (x k ,y k ). 

The solution of the BFKL equation obtained in the next section is also Mobius invari- 
ant in LLA and can be written in the form [37]: 



+°° r +oo 



[y 2 + n 2 /4) dv 



G vn {pi,p~2;pf,p^) 



[v 2 + (n - l) 2 /4] [v 2 + (n + l) 2 /4] uj - w(v, n) 



where for n — ±1 the integral in v is regularized as follows: 



+oo X 

dv — (p(v) = lim 

-no V 2 e^O 



- , 6(v 2 -e 2 ) 
dv 



<p(v)-2 



(114) 



(115) 



In the "energy propagator" {uj — u(v,n)) 1 the quantity u(v,n) is the eigen value of the 
BFKL equation [7]: 
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ufan) = ^ f 1 [x^l-^cos^ln x) - 1 

v ' ; 2vr 2 7o 1 - x L v ; 

L(i±M+i I /)-^(i)V (H6) 



2tt 2 V 2 
The Green function G un (p~i, p|; pj,p^>) is given below: 



Gvn(pi, p|; p$,p$)= J d 2 p E* n (p l/Q , p 2 > ) E m (p 10 , p 20 ) , (117) 



where 



^n(PlO, P20) = (0 | 0( Pl )0(p 2 )O ra (p O ) | 0) = (-^V)* , (H8) 

P10P20 P10P20 

are the solutions of the homogeneous BFKL equation [37]. They are equivalent to 
the Polyakov three-point function [25] for the case when the fields 0(pi) describing the 
reggeized gluons have vanishing conformal quantum numbers h and h. The composite 
field O vn describing the BFKL pomeron has the conformal weights 

' Tl — 1 Tl 

h= 2 +tu+ r h = 2 + lu ~2 (119) 

with real v and the integer conformal spin n in accordance with the fact that they belong 
to the basic series of the irreducible unitary representations of the conformal group. For 
g — > partial wave f w (114) is reduced to (110) because the functions E un (p 10 , p 2 o) have 
the following completeness property: 

_ _ +oc r+°° d v v 2 + n 2 /4 _ 
S 2 (pt -p?)5 2 (pl -p?) = Yl / —] [21 ^G un (pt,pt;pJ,pj). (120) 

n=-oo J -°° n IPll'l |P22'| 

For various physical applications of the four-point Green function (114) it is convenient 
to present G vn (p~i, pi; PiSP?) (H7) in terms of the hypergeometric functions taking into 
account its conformal invariance: 

G un (pl, p|; p7, p$) = ci x h x* h F(h, h, 2h; x) F(h, h, 2h; x*) 
+ c 2 X 1 ' h x* 1 ~ h F(l -h,l-h,2-2h; x) F(l-h,l-h,2- 2h; x*) , (121) 
where x is the complex anharmonic ratio: 



x = (122) 
P11' P22' 



and F(a, b, c; x) is defined by the series 



„, , \ ab a(a + 1) bib + 1) 9 

f(«Ac;*) = i + I ^* + V + . (123) 



The coefficients 
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c = K i ci _ _ r(2 - 2fe)r(2 - 2fe) [r(/^)r(/i)j 2 

° 2 2 tt 2 ' c 2 V* [r(l - h)F(l - h)f T(2h)T(2h) 
and the factor b 



(124) 



J n,v 



r(-iv + (i + 


n|)/2) r(ii/+ |n 


1/2) 


r(iv + (i + 


n|)/2)r(l -iu+ | 


n|/2) 



are obtained from ref.[37]. The ratio Ci/c 2 also can be fixed from the condition, that G 
is a single-valued function of its arguments. With the use of the various relations among 
the hypergeometric functions in Appendix we derive other representations for the Green 
function G, which gives us a possibility to continue it in the regions around points x — 1 
and x = oo . 

2.3 Operator product expansion of the ^-channel partial waves 

Let us consider now the properties of the solution (114) of the BFKL equation in the 
various interesting cases. It is convenient to use the mixed representation for the gluon- 
gluon scattering amplitude [37]: 



1 f i 

/*("J^> = J eX P(^ - ^ )~<f ) /"(Pi, P2', P~f, Pj) , P = Pl2 , P = Pi' 2 ' , 

(126) 

where R = ( Pl + p 2 )/2, R' = (p[ + p' 2 )/2. 

Because the dependence of the impact factors $ a ' b (pt, pi, if) from it according to 
(108) is simple: 

$ a < fe (p?, p|, It) = $ a (pH, if) exp(z^t) , (127) 
we have another representation for f w (q 2 ) (104): 

uq 2 ) = J -^f n ?, t) ait, 7) $ 6 (7; t) ■ (128) 

Here the function f%{~]?, p ) can be interpreted as the amplitude for the scattering of 

two composite objects with the sizes ~/f and p . Further, the impact factors $ a (~/7,lf) 

and $ 6 (p , if) describe the distributions over these sizes for fixed if. Note, that a simple 
space-time picture of the high energy scattering in QCD was developed in the framework 
of the dipole approach [40]. 

For small if the essential values of 7? and ~$ do not depend on if. As for ft*(~j?, p ), 
it has a weak singularity at small q 2 due to massless virtual gluons. Indeed, according to 

definition (122) for fixed ~/f, p and large R — R' ~ l/q the anharmonic ratio x is small: 

PP 

x- ~ -, /|2 < 1 



R-R 
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and therefore expression (121) for G vn can be simplified: 

G m (pt,pZ;p?,p2) ^ Cl x h x*~ h + c 2 x l - h x* l -~ h . (129) 
By putting it in eq.(114) we obtain, that at large distances R — R' the function / w equals 

PP 



fu(pi,P2;P?,P2) 

Here v{ui) is a solution of the equation 



(R-R') 2 



(130) 



u = u(u,0) (131) 

with Im{v) < 0. If analogously to (118) one will consider fu(p~i, p~2', frf, p^) as the four- 
point function of a two-dimensional theory: 

MpI, ti; pf, p?) = (o \(Kpl)(Kn) <Kp?)<Kp?) I o) (132) 

then its asymptotics at small p and p and fixed it! — R' is related with the anomalous 
dimension 

l(u) = ±(h + h) = ±+iv(u) (133) 
of the operators O w appearing in the operator-product expansion: 

0(pi)0(pl) ~ \p\2\ 1 O u (pt) , (f)(pf)(f)(pj) ~ |px' 2 ' | 27 O w (p?) • (134) 

The quantity 7(0;) is the anomalous dimension of the twist-2 operators constructed in a 
bilinear form from the gluon fields G tll/ (x). It can be calculated as an expansion in the 
parameter g 2 /u (cf. (26)): 

where ((x) = J2T=i k~ x i s the C-function. In LLA where 00 ~ g 2 the anomalous dimension 
is of the order of unity and it has a square-root singularity 



1= 2 + ^iuNM3) (136) 



at the point 



u = ^N c hi 2. (137) 

This value is the intercept of the Pomeranchuck singularity governing the asymptotics of 
the cross-sections a tot ~ s Wo . 

In the deep-inelastic regime of large Q 2 where pyy <C pn ~ pw one can again use the 
operator expansion for the product of the fields <f(pv) and <f(p2') and the asymptotics is 
governed by the anomalous dimensions 7(0;) (133) of the twist-2 operators with matrix 
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elements (118). As for the regime of large momentum transfers \t\ ^> |pi2p 2 ~ lpi'2'l" 2 ) 
we can use here the operator expansion of the product of the fields (p(pi) and (p(pi>), which 
corresponds to the asymptotic behaviour of f u (114) at x — > oo (see (122)). From the 
representations of G vn given in Appendix one can obtain the asymptotics of f u at large 
x: 

f u ~ c w In \x\ 2 , (138) 

where c u does not depend on p 12 and p 1/2 >. Therefore in the mixed representation we have 
at large \t\: 

lZ(j},T)^\c^ (139) 

which means, that in the momentum representation the function f w (k , k ) contains the 
contributions proportional to S 2 (ki) and 8 2 (k2). Using the normalization conditions for 
the wave functions of the initial and final particles, one can verify that at large momen- 
tum transfer the impulse approximation is valid and the hadron scattering amplitudes 
is expressed as a sum of the effective parton scattering amplitudes which do not contain 
infrared divergencies [39]. 



2.4 Asymptotic freedom and vacuum Regge poles 

Because the singularity of f u (p~i , pi; p~$ , py) at large it! — R' is weak there is a finite limit 
for the scattering amplitude at small q. Indeed, this amplitude can be written at the 
mixed representation in the form: 



I +oo f + oc 



y + n 2 /2) dv 



KniPi l)Evn(p', 0) 



[u 2 + (n - l) 2 /4] \v 2 + (n + l) 2 /4] u - n) 



(140) 

where the function E un (p, q) is the Fourier transformation of E un (pio, p 2 o) in po- In 
particular, for p q — > we obtain: 



Euo(p, q) ~ e(p) = ( |p| 2lv + 



q 2 p 



2iu 



exp i5(0, v 



(141) 



Here the phase 5(0, v) is the simple function of v (given in ref. [37]) which for v — > 
tends to the constant: 



5(0, i/) — > 7T . (142) 
The amplitude e(p) is a solution of the homogeneous BFKL equation at small p: 

2 N r 



^e(p) = axbP<9)e(p) , xiy) 

I 7T 



Re (V(1)-^(I + W )). (143) 



If we take into account the fact, that the QCD coupling constant is decreasing in the 
region of the large transverse momentum k ~ 1/p ^> q this equation should be modified 
by the substitution (see (13)) 
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The solution of the modified equation can be found easily. At not too small p where 
uj < ujq it can be written in a semi-classical approximation as follow [37] 



(7r lb r u (f^) , , \ 

4 + V{P ' U) ln ~ " •<» Xi) ) (145) 

where V(p, uj) satisfies the saddle point relation: 

^ln-^— = b X (v). (146) 
\p\ A z 

By matching the logarithmic derivatives of two above expressions (141) and (145) for e(p) 
at p ~ 1/q and using the approximate expression for x( v ) a ^ small v: 

Xiy) = — (2 In 2 - 7C(3)z/ 2 ) (147) 
we obtain the following quantum spectrum of V and a; at small v: 

, /3tt(A; + 3/4) ln2\ 1/3 . . 

c fe (-^ 2 ) = ^^ (2 1n2-7C(3)^(-^ 2 )) , (149) 
7T In -fa- 

where k — 0, 1, 2, .... Thus, as a consequence of the asymptotic freedom, the fixed square- 
root cut of the t-channel partial wave f% (140) situated at uj = uj q is substituted by an 
infinite sequence of the Regge poles. It is possible, that the residues of some poles including 
leading one with k — can be negligible small. For hadrons they are determined by the 
strong interactions at large distances. The dependence of the t-channel partial waves from 
the dynamics in the confinement region is even more significant for the fixed momentum 
transfer if 2 ~ A 2 . Because the intercept of the soft pomeron is small A ~ 0.08 <^ 1 [32] 
one can attempt to calculate its trajectory at small if in the perturbation theory [37]. In 
the region where uj a(l/ |p| 2 ) 1 the BFKL equation does not depend on uj and its 
solution has the simple form [37]: 



e(p) -cos (7 + 1/1 ln^— --] , a = b H dv> X W) (150) 
\4 \n A 2 uj J Jo 



vr , 1 

IpI 

where v 1 ~ 0.637 is the root of the equation x( v ) = an d a — 1-28 for rij = 3. It 
is resonable to expect, that one can neglect the ^-dependence of the equation also in 
the confinement region. In this case the small-p asymptotics of the soft pomeron wave 
function should have the form 

e(p) ~ cos ( 7 + 1/1 ln —i— - y?(lf 2 ) J • (151) 
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Here (p(lf 2 ) is a phase which is assumed to be a linear function (p(lf 2 ) = <p — ct with 
c > at small t = —if 2 . By matching above expressions for e(p) we obtain the Regge 
trajectories [37]: 

M-T)= ° p^ . (152) 
7r ft + </?( g ) 

Here without any restriction ip is assumed to be in the interval tt > <p > 0. To obtain 
positive values for the pomeron intercepts we should put k — 0, 1, ... and therefore the 
rightmost pole is situated at A > a/rr = 0.4 [37]. With increasing t = — (if) 2 the function 
(f(lf 2 ) is decreasing and can reach the point ip — where according to (152) the first 
Regge pole goes to the infinity. Near this point one should use more accurate expressions 
leading to a more moderate growth of uo(t). Nevertheless, to avoid the non-physical 
behaviour of the scattering amplitude at positive t the residue of the first Regge pole 
should be small for all momentum transfers. In an analogous way the second pole with 
k = 1 also could have a small residue. For the position of the next poles the value of 
the unknown phase tp is not essential. For example, the intercept of the third Regge pole 
approximately equals to the phenomenological intercept of the soft pomeron. 

Note, that in the momentum representation at t — the solution of the BFKL equation 
with a fixed QCD coupling constant can be written in the form: 

f(^k,k') = 7 ±- 2 / du £ {k/k) ? \ k/k } , (153) 



2tt 2 J-oo n fr' 00 u-u(v,n) 

where k and k! are the transverse momenta of incoming and outgoing particles at the 
t-channel. This amplitude is normalized in such way, that at small g we obtain: 

f(u,k,k') -> - 5 2 (t -t) 

Near the leading singularity u = u the region of small v is essential and we obtain the 
well known diffusion expression for the virtual gluon cross-section at high energies [7] 

/ ,9 ,/9n r du ,, , ... s W0 . ln 2 (|A;| 2 / l^| 2 )s , , 

a(s, k , k ) ~ / — f(u, k, k ) ~ ~ 7 == exp 11 7 1 U 154 

J 2tti Vain s c In s 

where 

c = 56^C(3) 

71 

is the diffusion constant. 



3 Multi-Regge processes in QCD 

For the deep-inelastic scattering at small Bjorken variable x the gluon distribution g(x, k±) 
depending on the longitudinal Sudakov component x of the gluon momentum k and on 
its transverse projection k± in the infinite momentum frame of the proton \~J?a\ ~^ 00 can 
be expressed in terms of the imaginary part of the gluon scattering amplitude at t — 
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in the Regge regime of high energies y/s = y/2pxp b an d fixed momentum transfers. The 
most probable process at large s is the gluon production in the multi-Regge kinematics 
for final state particle momenta k = pa>, k\ = q\ — q 2 , ■■■k n = q n — q n +i, k n+ \ = pb> (see 
Introduction): 



i— 1 n+1 n ^ 2 ^ 2 

s > Si = 2ki-iki > U = q] = (pa~Y1 ^) 2 , II s * = s II ^ , k]_ — -k . (155) 

r=0 i=l i=l 

In LLA the production amplitude in this kinematics has the multi-Regge form [7]: 

n+l 

A L 2 t A 2+n = A\Z\ +n II ^ ■ (156) 
i=i 

Here s^ 1 ^ are the Regge-factors appearing from the radiative corrections to the Born 
production amplitude A 2 r ^ 2+n . The gluon Regge trajectory j — 1 + cj(t) is expressed in 
terms of the quantity: 

<"(*) = "0 / ^2,^^ , t = -If 2 • (157) 



fc (Y- k) 2 



Infrared divergencies in the Regge factors cancel in a tot with analogous divergencies 
in the contributions of real gluons. The production amplitude in the tree approximation 
has the following factorized form [7] 

A 2 Z e 2+n = 2^r 1 i^ri 1 i....^ 1 ^r^ 1>B ^^ 1 r 2 . (158) 

II f2 In+1 

Here A, and A', B', d r . (r = 1, 2...n) are colour indices for initial and final gluons corre- 
spondingly. T£ b = —if a bc are generators of the gauge group SU (N c ) and g is the Yang-Mills 
coupling constant. Further, 

Ti = \ey;*T^\ T r r+hr = ~r,(q r+1 ,q r )e^(k r ) (159) 

are the reggeon-particle-particle (RPP) and reggeon-reggeon-particle (RRP) vertices cor- 
respondingly. The quantities A r = ±1 are the s -channel gluon helicities in the cm. 
system. They are conserved for each of two colliding particles: T 1 = 5y\, which is not 
valid in the one loop approximation [41]. The tensor Y uu ' can be written as the sum of 
two terms: 

= _ ^( n +y]_( n +Y, (160) 
Pa 

where we introduced the light cone vectors 

n- = ^,n + = E = Vs/2, n+n' = 2 (161) 
E E 

and the light cone projections = k a n^ of the Lorentz vectors k a . The first term is the 
light cone component of the Yang-Mills vertex: 

l uu ' + = {p + a+P + a>W u ' - 2pA(n + y - 2p» A ,(n + ). (162) 
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The second (induced) term in (160) is a coherent contribution of the Feynman diagrams 
in which the pole in the t-channel is absent. Indeed, it is proportional to the factor q 2 
cancelling the neighbouring propagator. 

Similarly the effective RRP vertex T(q 2 , qi) can be presented as follows [7] 

T°(q 2 , q x ) = - + 2g 2 2 ^, (163) 

where 

7 CT +- = 2q 2 T + 2q1 - 2(n~) a kf + 2(n + ) a k^ (164) 

is the light-cone component of the Yang-Mills vertex. 

Due to the gluon reggeization the above expression for the production amplitude in 
LLA has the important property of the two-particle unitarity in each of the ti channels 
(see (74)). 

Furthermore, it satisfies approximately the unitarity conditions in the direct channels 
Si with the intermediate particles being in the multi-Regge kinematics [7, 37]. Note, 
that from the general analytic properties of the inelastic amplitudes in the multi-Regge 
kinematics discussed in Introduction (see (86, 87)) one can obtain dispersion relations 
in a differential form. Indeed, in accordance with the fact, that in LLA the real part of 
the production amplitudes for the negative signature is much bigger than their imaginary 
part we can simplify the signature factors for different dispersion contributions near the 
point j — 1. Using this simplification one can derive the following relation 

d 

Im Sai A 2 ^2+n = tt^t-t — -ReA 2 ^ 2+n . (165) 
j om{s a i) 

Here Im Sa .A 2 ^ 2+n means the imaginary part of the production amplitude in the channel 
s ai = ipA' + Dr=i K) 2 - Each of these imaginary parts can be calculated through the prod- 
ucts of Re A 2 ^ 2+ i summed over all intermediate states and integrated over the produced 
particle momenta with the use of the unitarity conditions. Finally one can verify that 
the above multi-Regge production amplitudes in LLA satisfy the obtained "bootstrap" 
equations [7, 37]. 

Note, that T a has the important property: 

(hrr^q 2 , qi ) =0, k 1 = q 1 -q 2 , (166) 

which gives us a possibility to chose an arbitrary gauge for each of the produced gluons. In 
the left (/) light cone gauge where pAe l (k) = the polarization vector e l (k) is parametrized 
in terms of the two-dimensional vector e l ± 

e l = e\ - k -^p A (167) 
kpA 

and satisfies the Lorentz condition k e l = 0. The matrix element of the reggeon-reggeon- 
particle vertex T takes an especially simple form [42] 

rl tl = Ce* + C*e, C=4?> ( 168 ) 

K 1 
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if we introduce the complex components 

6 — — &x ~\~ 1&y •, 6 — — ) ^ ~~ ^x ~~ l~~ ^fcy 5 k — — kx iky (^169) 

for transverse vectors e±, k±. The factors g* and q<i in expression (168) guarantee the 
vanishing of the inelastic amplitude at small momentum transfers, which is a consequence 
of the fact that the virtual gluons interact in this limit with the total colour charge, whose 
matrix elements are zero between the states with the different number of gluons. The 
singularity 1/k* in C reproduces correctly the bremsstrahlung factor for the production 
amplitude in the soft gluon emission theorem. 

The above complex representation was used in [42, 43] to construct the effective scalar 
field theory for multi-Regge processes which in particular can be applied for the derivation 
of the equations for compound states of several reggeized gluons [44, 45]. The correspond- 
ing effective action was obtained later from Yang-Mills one by integrating over the fields 
describing the highly virtual particles [43]. 

The effective action describing multi-Regge processes can be written in the form in- 
variant under the abelian gauge transformations 8V£ = id^x a f° r the physical gluon fields 
Vfj, provided that the fields A± corresponding to the reggeized gluons are gauge invariant 
(SA : -0): 

S mR = J d A x{\{F; v f + \{d ±a AX){d ±(T A«_)+ 
l - g [-A% {F_ c T a idZ l F_ c ) - A a _(F +a T a id' 1 F +a ) + {dZ 1 F» a ){A_T a id a A + ) + 

+ (d?F%)(A + T a id A-) + t(j-^Fl_)(d a A + )T a (d a A^) + iF°_(A_T a A + )}} , (170) 

where F^ v = d^V v — d u V^ and we introduced the light-cone components of the Lorentz ten- 
sors in accordance with (161). The fields A± satisfy the kinematical constraints d±A T = 
equivalent to the condition that in the multi-Regge kinematics the reggeized gluon takes 
the negligible part of energy from the colliding particles. The Feynman vertices which are 
generated by this action coincide on the mass shell with the effective vertices (160, 163) 
for the reggeon-gluon interactions. However, for virtual gluons these vertices are different, 
which leads to the inconsistency of this theory. Such drawback is absent in a more general 
nonabelian effective field theory which will be discussed in the next section. 

3.1 BFKL Pomeron 

Using the explicit expressions for production amplitudes in the multi-Regge kinematics 
one can calculate the imaginary part of the elastic scattering amplitude in LLA with the 
pomeron quantum numbers in the crossing channel related with the high energy asymp- 
totics of the total cross-section. Its real part is small in accordance with the fact, that this 
amplitude has the positive signature. Due to the factorized form of the production ampli- 
tudes one can write down the Bethe-Salpeter equation for the vacuum t-channel partial 
wave f u describing the pomeron as a compound state of two reggeized gluons [7]. The 
convenience of the (^-representation is related with the angular momentum conservation 
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in the t-channel. The contribution to the integral kernel of the BFKL equation from the 
real gluons is proportional to the product of the effective vertices calculated in the light 
cone gauge [42]: 

C(p 1 ,pv)C*(p 2 ,p 2/ ) + h.c. = PlP ^ 2 ' + h.c. (171) 

where pi,p 2 and pr,p 2 ' are the corresponding complex transverse components of initial 
and final momenta in the t -channel (q = pi + p 2 = py + p 2 >). In turn, the contribution 
related with virtual corrections to the production amplitudes is proportional to the sum 
of the Regge trajectories of two gluons: 

u(-pt 2 ) + uj(-p?) ~ In | Pl | 2 + In | p 2 I 2 +c, (172) 

where the constant c contains the infrared divergent terms which are cancelled with the 
analogous terms from the real contribution after its integration in k. The final homoge- 
neous equation for the wave functions of pomerons being the singularities of f u takes the 
form [46] 

= H 12 V, E = ~^r. (173) 

g 2 N c 

Here the " Hamiltonian" H 12 is [7, 47] 



H 12 = In | Pl | 2 + In | p 2 | 2 + 1 -{p\p 2 In | p 12 | 2 p lP * 2 + h.c.) - 4^(1) (174) 

I Pl | | P2 \ 

where i/j(x) = ^lnT(a;) and T(x) is the Euler T-function. In the above expression 1/pt 2 
are the gluon propagators. We introduced the complex components pu = %k + Wk for the 
impact parameters canonically conjugated to the momenta pk = i-^ and performed the 
Fourier transformation: 



A: 



12 



In | P12 I 2 , (175) 



where pik = p% — Pk- The expressions 



In | Pi | , \Pi\~ 1 (176) 

are the integral operators in the impact parameter representation. The Hamiltonian (174) 
has the property of the holomorphic separability [46]: 

H 12 = h 12 + h* 12 , E = e + e, (177) 

where e and e are the energies correspondingly in the holomorphic and anti-holomorphic 
subspaces: 

e^(pi,p 2 ) = h 12 i/j(p!,p 2 ) , eTp(pl,p 2 ) = h* 12 $(pl,p* 2 ), , p|) = i/j (178) 
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The holomorphic hamiltonian is [47] 

h 12 = — In (p 12 ) Pl + — In (p 12 ) p 2 + ln(p lP2 ) - 2^(1). (179) 

Pi P2 

One can verify the validity of another representation for H 12 : 

hu = P12 ln(pip 2 ) Pu 1 + 2 ln(p 12 ) - 2^(1). (180) 
Further, using the following identities: 

2 In d + 2 In p = i^(-pd) + ij)(dp) = ip(-pd) + pd) (181) 

and 

2 ln(p 2 <9) - 2 ln(p) = ^(p<9) + ^(-p 2 dp~ l ) = ^j(pd) + ^(1 - pd) , (182) 
we can derive the following formulas [47] : 

h = Hp 2 uPl ) + \n(p 2 12 p 2 ) - 21n(p 12 ) - 2^(1), (183) 



h = ^(Pndi) + ^(P2id 2 ) + ^(1 + pai^i) + ^(1 + Pnd 2 ) - 2^(1) . (184) 

Here V'Oe) = r'(x)/r(x) and T(x) is the Euler T-function. From eq. (183) one can easily 
verify that h is invariant under the Mobius transformations: 

aPk + b 

Pk -> . j 185 

cpfe + rf 

for arbitrary complex values of a, 6, c, d. It means, that solutions of the homogeneous 
BFKL equation belong to irreducible unitary representations of the Mobius group. The 
generators of this group for a general case of n-particle interactions are 

Tl TL Yl 

M z = J2 Pk d k ,M~ = J2d k ,M + = -J2p 2 k d k . (186) 

k=l k=l k=l 

Its Casimir operator is 

M 2 = (M z ) 2 - -(M M ■ M M ) = -Y,Pr S drd s . (187) 

For the wave function of two particles in the holomorphic subspace we can use the Polyakov 
ansatz (cf. (118) with the substitution h — > m): 

<Mpio, P20) =< | y?(piMp 2 )O m (p ) | >= {-E^-) m . (188) 

P10P20 

The conformal weight m = \ + iv + ^ ot the composite operator O m is related with its 
anomalous dimension d — 1/2 + iv and its conformal spin n. This operator belongs to 
the basic series of the unitary representations provided that v is real and n is integer. In 
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turn, the fields <p(pi) describe the reggeized gluons and have the trivial quantum numbers 
d = n = 0. The holomorphic factor ip m is an eigen function of the corresponding Casimir 
operator: 

M 2 ^ m = m{m - (189) 
Simultaneously it is an eigen function of the BFKL equation in the holomorphic subspace: 

h^ m = e^j m , e = ^(m)+V(l-m)-2^(l). (190) 

The eigen value e can be obtained from representation (184) for h if one would integrate 
the both sides of the corresponding Schrodinger equation in (178) over the coordinate p 
with the use of the relation: 

/ dp ^(p 10 ,p 20 )~ pl~ m . (191) 

The second Casimir operator M 2 * is expressed through the conformal weight m — \ + 
iv — |. The total energy is 

E = ip(m) + -m) + ip(m) + - m) - 4^(1). (192) 
After simple transformations one can rewrite E as follows 

E = 4Reip{- + -4^(1). (193) 

The minimum of the energy is obtained for v = n = and equals Eq = —8 In 2. There- 
fore the total cross-section calculated in LLA grows as s (where A= (g 2 N c /n 2 ) In 2), 
which violates the Froissart bound a tot < cln 2 s [7]. One of the possible ways to improve 
LLA is to use the above effective field theory for multi-Regge processes [42, 43]. 



3.2 Compound states of reggeized gluons in multi-colour QCD 

The simple method to unitarize the scattering amplitudes obtained in LLA is related with 
the solution of the BKP equation [44] for compound states of n reggeized gluons: 

Em = Y J H lk ^>. (194) 

i<k 

Its eigen value E is proportional to the position uj = j — 1 of the singularity of the t-channel 
partial wave: 

8tt 2 

e = -jw"- (195) 

Note, that a non-trivial example of the BKP dynamics is the equation for the Odderon 
which is a compound state of three reggeized gluons [36, 45]. 

The pair Hamiltonians H ik in (194) have the property of the holomorphic separability: 

H tk = -^(h tk + h* k ). (196) 
33 



The group generators act on colour indices of the gluon i. The holomorphic pair 
hamiltonian is 

1 1 

h ik = — \n(p ik )Pi H \n(p ik )p k + ln(pip k ) - 2^(1) . (197) 

Pi Pk 

Similar to the Pomeron case we introduce the complex coordinates p k = x k + iy k (k = 
1,2, ...n) and their canonically conjugated momenta p k = i g^j in the impact parameter 
space (note, that p ik = pi — p k ). The above Schrodinger equation for \1> is invariant [46] 
under the Mobius transformations: 

/ a p k + b 

Pk = ; 

c p k + a 

for any complex values of a, b, c, d. 

According to t'Hooft for the multi-colour QCD (N c — > oo) only planar diagrams in the 
colour space are important. Because the colour structure of the eigen function at large 
N c is unique, due to eq. (196) the total hamiltonian H can be written as a sum of the 
mutually commuting holomorphic and anti-holomorphic operators [46]: 

H = hh + h*), [h,h*] =0. (198) 

The colour factor 1/2 appears because at large iV c the neighbouring gluons are in the 
octet state. The holomorphic hamiltonian is 

n 

h = '£h i , i+1 . (199) 

i=i 

Thus, in the multi-colour QCD the solution of the Schrodinger equation (194) has the 
property of the holomorphic factorization: 

* = E C kMpl, -Pn) Mpl, -P* n )- (200) 

where ip and ip are correspondingly the analytic and anti-analytic functions of their argu- 
ments and the sum is performed over all degenerate solutions of the Schrodinger equations 
in the holomorphic and anti-holomorphic subspaces: 

eip = hip, e* ip* = h* ip*, E = ^(e + e*). (201) 

These equations have the nontrivial integrals of motion [47]: 

t(9)=trT(9), [t(u),t(v)} = [t(9),h] = 0, (202) 

where 9 is the spectral parameter of the transfer matrix t(9). The monodromy matrix 
T{9) is constructed from the product of the L-operators 

T(9) = L 1 (9)L 2 (9)...L n (9) (203) 

expressed in terms of the Mobius group generators: 
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L k (6)=( 9 + ip 2 f R ldk A. (204) 
v ' \ -ip%d k 0-ipA J v ' 

The solution of the Schrodinger equation (201) in the holomorphic subspace is reduced to 
a pure algebraic problem of finding the representation of the Yang-Baxter commutation 
relation [47]: 

T hi[ (u)T i2i/2 {v){v -u + iP 12 ) = {v-u + iP 12 )T i2i , 2 (v)T hi[ { U ) , (205) 

where the operator P\ 2 in its left and right hand sides transmutes correspondingly the 
right and the left indices of the matrices T(u) and T(v ). Moreover [48], Hamiltonian (199) 
coincides with the local Hamiltonian for a completely integrable Heisenberg model with 
the spins belonging to an infinite dimensional representation of the non-compact Mobius 
group and all physical quantities can be expressed in terms of the Baxter function Q(X) 
satisfying the equation: 

t(X)Q(X) = (A + i) n Q(\ + + — i) n Q(\ - i), (206) 

where t(\) is an eigen value of the transfer matrix. The solution of the Baxter equation 
is known for n = 2. In a general case n > 2 one can present it as a linear combination of 
the solutions for n = 2 with a recurrence relation for the coefficients d k . For n = 3 this 
relation takes the form: 

Mk(A) = 2{2k + l) {k - m + 1 )( k + rn)(d k+1 (A) + d k ^(A)) (207) 

with the initial conditions do — 0, d± — 1. If one will consider for simplicity the inte- 
ger values of the conformal weight m the quantization condition for eigen values A is 
d m -i(A) = 0. Although the orthogonality and completeness conditions for the polynomi- 
als d k (A) are known: 

m ~ 2 2(2k + 1) 

E Kk + m \ m+ \ ){k + m) MA) - 5 aA _ i{ A) , A* ; (208) 

d k (A)d~ k (A) fc(fc + l)(fc-m + l)(fc + m) d 

§ - 6kk mTT) ' dk{A) - dA dk{A) ' (209) 

their full theory is not constructed yet. It does not allow us to calculate analytically the 
intercept and wave function of the Odderon in QCD [45, 47]. 



4 Effective action for small x physics in QCD 

All above results are based on calculations of effective Reggeon vertices and the gluon 
Regge trajectory in the first nontrivial order of perturbation theory. Up to now we do not 
know the region of applicability of LLA including the intervals of energies and momentum 
transfers fixing the scale for the QCD coupling constant. One should develop also a self- 
consistent approach for the unitarization of the BFKL pomeron. Therefore it is needed 
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to generalize the effective field theory of ref. [42] to the quasi-multi-Regge processes in 
which the final state particles are separated in several groups consisting of an arbitrary 
number of gluons and quarks with a fixed invariant mass; each group is produced with 
respect to others in the multi-Regge kinematics. The production of two particles with a 
fixed invariant mass is the simplest example of such processes [49]. 

At high energies the rapidity y = \ In |p constructed from the light-cone compo- 
nents /c ± = k a n^ of the particle momenta is similar to the time in quantum mechanics. 
The corresponding Hamiltonian is determined by the interaction of gluons and quarks 
with a nearly equal rapidity. Let us introduce the parameter 77 much smaller than Ins. 
The gauge-invariant effective action S e ff local in a rapidity interval (yo — r],yo + rj) was 
constructed recently [51] and includes apart from the usual Yang-Mills action also the 
reggeon-particle interactions (cf. (170)): 



S eff (v,A ± ) = -j d'xtr [\Gl(v) + (A_(v) - A_)j^ + (A + (v) - A + )f_ e ° 



(210) 



where the anti-hermitian SU(N C ) matrices v a and A± describe correspondingly the usual 
and reggeized gluons. Because the action is local in the rapidity space, we omit temporally 
y as an additional argument of these fields. The reggeon current j± 9 is expressed in terms 
of A± as follows: 

J± 9 = dlA ± , (211) 

which guarantees, that the interaction disappears on the mass shell k 2 = 0. 
The fields A± are invariant 

5A ± = (212) 
under the infinitesimal gauge transformation 

5v a = [D a , X ] (213) 

with the gauge parameter x decreasing at x — > 00, but they belong to the adjoint repre- 
sentation of the global SU (N c ) group and are transformed at constant x as follows: 

5A ± =g[A ± , X ] • (214) 



As usually, 



Gp,(v) = ^ [Dp, D u ] = d^v u - dyV^ + g v v \ , D a = d a + gv a . (215) 



The fields A± obey the additional kinematical constraints 

d + A_ = , d-A + = (216) 

meaning that the Sudakov components a, (5 of the reggeon momentum q = aps+P PA+k± 
are negligible small in comparison with the corresponding big components fik-i of 
the neighbouring particle momenta. Such simplification takes place at the quasi-multi- 
Regge kinematics where the gluons in the final and intermediate states are separated 
in several clusters. The invariant mass of each cluster is restricted from above by a 
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value proportional to exp(rj) and the neighbouring clusters are significantly different in 
their rapidities: yu-i — Uk ^> V- Further, the Sudakov components (5k of their total 
momenta are strongly ordered: at ^> cxk-i, Pk Pk-i and the transverse momenta 
k± are restricted. The effective action describes the self-interaction of real and virtual 
particles inside each cluster and their coupling with neighbouring reggeized gluons. Note, 
that because small fractions of the Sudakov parameters are transmitted from one to other 
clusters the constraints d±A T = (216) are not absolute. The composite reggeon field 
A±(v) is given below 

A±(v) =v±- gv±-^-v± + g 2 v±-^-v±-^-v± - ... (217) 
o± o± o± 

and can be written in the explicit form 

A ± (v) = v^di = --d ± U(v ± ) , U(v ± ) = Pexp(-§ f dx'^x')) , (218) 

g z J -co 

where the integral operator D± 1 d± is implied to act on an unit constant matrix from the 
left hand side and the symbol P means the ordering of the fields v in the matrix product in 
accordance with the increasing of their arguments x !± . We chose for the Wilson exponent 
U(v±) the following boundary condition: U — > 1 at x — > — oo. This exponent does 
not have infrared singularities in the essential integration region in accordance with the 
fact, that the operators dg 1 are restricted from above for all particles inside the cluster. 
Because j reg (211) in the momentum representation contains the factor t = q 2 killing the 
pole in the neighbouring reggeon propagator, the corresponding scattering amplitudes 
satisfy the Steinman relations forbidding simultaneous singularities in the overlapping 
channels t and s. The interaction terms describe contributions of the Feynman diagrams, 
in which the gluons in the given rapidity interval (y ~ Vi Ho + rj) are coherently emitted 
by the neighbouring particles with essentially different rapidities. 

The interaction terms of the action are gauge invariant due to the following relations 

At, fi d ] = , 3 l " d = ^ tr ( A ±(v) d 2 a A T ) (219) 

where j™ d is the induced current. For the particles belonging to the same cluster the 
parameter r] plays a role of the ultraviolet cut-off in their relative rapidity. 

Note, that action (210) is similar to the Legendre transformation for which the A±- 
dependence of the current j± g would be fixed by the stationary condition 5 F = for the 
free energy: 

expHF(j; e9 ,^±)) = / Udv(x) exp (-iS eff (v,A ± )) (220) 

with taking into account the Faddeev-Popov ghosts. Below some arguments supporting 
our choice (211) of j± 9 will be given. Because we do not use the Legendre transformation, 
there is a finite renormalization of the kinetic term for the fields A± in (210) even in the 
tree approximation: 

- A-ft* - A +3 L ea - (dtjA— ) (d a A + ) (221) 
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related with the transition A — > v — > A induced by the bilinear terms in the gluon and 
reggeon fields. 

4.1 77 - independence 

The physical results should not depend on the auxiliary parameter 77 due to cancellations 
between the integrals over the invariant masses of the produced clusters and the integrals 
over their relative rapidities (for which 77 plays a part of the infrared cut-off). This criterion 
is very important for the self-consistency of the effective action (210). For big values of 
77 the most essential polynomial contribution from integrals over particle rapidities inside 
each cluster corresponds to the production of an arbitrary number n of subgroups of 
real and virtual particles in the quasi-multi-Regge kinematics. We show below, that 
the interaction of the particles inside the whole cluster can be described by the same 
effective action for each subgroup with a smaller cut-off and the integrals over their 
relative rapidities reproduce correctly the contribution proportional to the polynomial 

Pn-M- 

To begin with, let us consider only two subgroups of real and virtual particles with 
the fixed invariant masses being restricted from above by a value proportional to exp(7? ) 
where 770 is a new parameter which is assumed to be much smaller than initial one 77 
for the whole cluster. The rapidities y l and y 2 of these subgroups satisfy the condition 
<y l — y 2 < 2i] leading to the integral contribution proportional to 77. 

The field v a can be presented as the sum of three Fourier components: 

v = v l + v 2 + v 12 . (222) 

The fields v 1 and v 2 describe the particles inside two extracted subgroups and v 12 is 
related with the strongly virtual gluons having the relatively big Sudakov parameters 
(3 ~ k + and a ~ k~ of the order of value of (3\ ~ k± and a 2 ~ &2 f° r the first and second 
subgroups correspondingly. Because on and /3j are strongly ordered in the quasi-multi- 
Regge kinematics: «2 > «i, A > P2, the virtuality of fields v 12 is large: M 2 = Saf3 S> k\. 

Let us consider the functional integral over the heavy fields v 12 for fixed v 1 and v 2 
(cf. [43]). In expression (217) the integral operators dg 1 correspond to the propagators of 
the strongly virtual neighbouring particles producing gluons within the chosen rapidity 
interval (yo—r], yo+rj). Because these virtualities are significantly different for the emission 
of gluons inside above two subgroups we can simplify the composite fields A±(v) by 
neglecting the more virtual fields: 

AJv) ~ --U (v 2 _) + ... , A + (v) ~ -?± U(vV) + (223) 

9 9 

The factors U(v\) and U(v^_) in the above expansion lead to the corresponding inter- 
action terms A + (v 1 ) d 2 A- and A^.(v 2 ) d^,A + in the effective Lagrangians L 1 and L 2 for 
each of the above two subgroups. To get the other two terms in these Lagrangians one 
should calculate the functional integral over v 12 . It leads to the contributions depending 
only on v\ and v\ because they enter with the coefficients proportional to the big Su- 
dakov parameters a<i ~ k£ or f3\ ~ /cf . Due to the gauge invariance of the initial action 
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these contributions should be also gauge invariant. This property fixes their form almost 
completely: AL 1 = A-(v}_) d 2 a + and AL 2 = A + {v\) d\a_ up to two arbitrary functions 
a±(x) satisfying the kinematical constraints <9_a+ = <9 + a_ = 0. Further, one can integrate 
the obtained expression exp(i5'^ n * + iS™ t ) over a± with an arbitrary weight depending on 
a± but to achieve an agreement with the perturbation theory the total Lagrangian should 
contain at least the quadratic term — 2 (d a a + ) (<9 a a_) leading to a correctly normalized 
propagator of fields a± responsible for the interaction of the particles with the relative 
rapidities smaller than 77 (and bigger than the new intermediate parameter 770). 

The integration over a± with the subsequent integration over the relative rapidity gives 
the result proportional to rj in accordance with our expectations. Non-linear interactions 
of the fields a± would lead to the contributions having higher powers of rj and therefore 
they should be absent for the two-group kinematics. Note, that these higher power con- 
tributions appear if we shall consider several subgroups of real and virtual particles in the 
quasi-multi-Regge kinematics inside the given cluster. In this more general case one can 
also verify, that the integration over momenta of the strongly virtual particles gives the 
interaction terms linear in a± for each subgroup. It means, that the functions a± coincide 
in fact with the fields A± and describe the interaction between the subgroups of particles 
within the additional interval (770, rj) of their relative rapidity. Thus, in the framework of 
the effective action approach the 77 -dependence of integrals over particle momenta inside 
each group is compensated by the 77-dependence of integrals over the relative rapidities of 
these groups. In particular for 77 = y the fields A± are absent and we return to the initial 
Yang-Mills theory. 



4.2 Classical equations 

One can verify, that the composite fields A±(v) are changed under the infinitesimal trans- 
formations of the gauge group 5v a = [D a , x] according to the abelian law: 5A±(v) = d±ip. 
Therefore after integrating by parts and using the kinematical constraint d±A T = we 
conclude, that the effective action S e ff is gauge invariant for \ decreasing at 00 provided 
that the reggeon fields are fixed: 8A± = 0. For constant x the action is also invariant 
because the fields A± belong to the adjoint representation of the gauge group SU(N C ). 

The effective action S e ff has a nontrivial stationary point v — v satisfying the following 
Euler-Lagrange equations [51]: 

= 0, (224) 
where the induced current f™ d equals 

ft = j^d T (d 2 ±a A ± ) d T ^- , jZ d = (225) 

and due to (219) satisfies the covariant conservation law: 

= . (226) 

We remind, that the integral operators D^- l d T and d^D^ 1 are implied to act on unit 
constant matrices correspondingly from the left and right hand sides (in the second case 



[D, 



At, G, 



up 



D 7" 
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after integrating by parts the operators d T should be substituted by — d T acting on A±) 
Using the following expression for the Yang-Mills current 



J™ = - [A,, = 

d^daVa - dfcvp + g([v a , d„v a ] - [d a v a , u„] -2[v a , d^v^]) - g 2 [v a , [v a , v^]] 
we can rewrite the above equations in the form: 

where 



3? = -9 



v a ,G, 



and 



j± = -g [ v cr, G a ±] + j± d . 

The above equations are self-consistent because the current j is conserved: 



ind 



+ 



n Air. 



ind 



0. 



(227) 

(228) 

(229) 
(230) 

(231) 



We can construct a perturbative solution v — v of the classical equations. For example 
in the Landau gauge where 

d a v a = 



one obtains 



and 



^ = EsX, f± = A±, fL = o. 



n=0 



(232) 
(233) 

(234) 



?1=0 



The higher order coefficients /" , A™ satisfy the reccurence relations: 



n-l 



3 2 f n = A n - V 

^aJ a a / t 



i=0 



i-1 



3=0 



and 



n-l 

= - e > A 

i=0 



(235) 



(236) 



where we used the gauge properties (219) of j ± d . 

The coefficients A" can be presented in the explicit form: 

n i n 

^l = J2(— d ^9lA ± (d T —) n . k , A™ ± = 



where 



1 OO 1 oo 



T perm r =0 r=0 



(237) 



(238) 
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1 OO 1 oo 



T permr=0 r=0 



(239) 



In these expressions the summation is performed over all permutations of the non-commuting 
integral operators (— ^-/^) and(— f^-^) correspondingly and the following reccurence re- 
lations are valid: 



1 n_i 1 

M 7r a F )n = -Eir 1 " i (7r^ i ' 

T i=0 T 

Maf^-)n = -E(a,=^-)i^- 1 - i - 



(240) 
(241) 



Using equations (235, 236), we can construct / ra and A n in several orders of pertur- 
bation theory taking into account that dzpA± = 0. Namely, 



f°± = A ± , / CT 0± = 0, A° ± = %A ± , A^ = 0; 
= - I [A T , d ± A ± \ , Ai = \d 2 a A ± , d- l A T 



+ 2 



A-,d^A, 



and 





+ 




+ 


fl,d a AJ 


+ 





A + ,drt - \A.,d + fi\+n a , 



where fi^r = and 



(242) 
(243) 
(244) 



(245) 



ih - A* - 2 fld a A ± ] - - [At, [A T , A±\] , (246) 

(247) 

These perturbative terms /" satisfy the Landau gauge condition d a f™ = 0. For the quasi- 
elastic kinematics when A + = (A- = 0) the solution of the reccurence equations is 
trivial: V- = A^ (v+ = A + ). 

Note, that the Euler-Lagrange equations for the effective action S e ff can be obtained 
also from the usual Yang-Mills action in the form 



A± 



d-Ulf'dlAJ + U^Al 



[D T ,[D (T ,G (7± ]]=0, [D a ,Gj p ] = 



(248) 



if one will restrict oneself to the quasi-gauge variations 8v± = [D±,x±] under which the 
interaction terms are stationary. Therefore the reggeon fields A± can be considered as 
some parameters of the solutions of eqs (248) and the effective theory can be derived 
from the gluodynamics by imposing certain constraints for the initial Yang-Mills fields: 
A±(v) = <p± where ip± is an abelian field. 
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4.3 Reggeon action in the semi-classical approximation 

To obtain the reggeon action in the tree approximation we should substitute v c — > v c in 
S e ff [51]. The interaction terms in S e ff can be expressed in terms of the Wilson integrals 
with the contours displaced along two light-cone lines: 

Slf f = i J d 2 x ± tr (^J°° dx + {T{v^) - l)d 2 a A + + J°° dx~{T{v + ) - 1) d 2 a A^j , 

Q f + OO I 

T(v ± ) = P exp (-§ / dx±v ± ) = lim — d ± , (249) 

where we used the constraints d±A T = 0. Note, however, that because this constraint is 
not absolute and T(v±) contains singular operators d± l there could be uncertainties. In 
this case the initial representation for the action is preferable. 

The expression d±Dg 1 appearing in expression (225) has the following integral repre- 
sentation: ^ 

d±Dl 1 = P exp (-- / dx' ± v±). (250) 
2 Jx ± 

Therefore it tends to 1 at x ± oo and the interaction term can be expressed in terms of 
the induced current j± d (225): 

S™ = - f d 2 x ± dx + tr (f+ d (oo) - d 2 a A + ) + - f d 2 x ± dx~ tr (j md (oo) - d 2 a A^) . (251) 
9 J 9 J 

For this choice of the boundary conditions for the perturbative coefficients A!^ of j± d the 
factors (-^-d T ) k and (d T -^-) k are defined as the coefficients of the perturbative expansion 
of the following expressions: 



1 oo / n+1 fX T \ 



(252) 



1 oo / n +l „oo \ 

d T — = P J] exp (-^- da*%) . (253) 

In the particular case of the quasi-elastic process, when there are only two clusters, 
after integration over the fields A± we obtain for the effective action 

STff = ~\j d^xtrGlM +AS, (254) 
where A S is the action for a two-dimensional a-model: 

A 5 = -j 2 J d 2 x ± tr (d ±a T(v_)) (d ±a T(v + )) , (255) 

introduced firstly by E. Verlinde and H. Verlinde using other arguments [52]. In our 
approach the Yang-Mills interaction of the fields v inside of two separated clusters is also 
essential. 
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By inserting the perturbative solution v a = v a , j± = j± of classical equations in 
S e ff one can write the reggeon action in the tree approximation as follows: 

5*7/ = - J d A x tr (s 2 + gs 3 + g 2 s 4 + ...), (256) 

where 

s 2 = d^A + d^A_ , (257) 

s 3 = ~(d 2 a A-) A + d- 1 A + - (dlA + ) A_dZ l A_ , (258) 

s± = -{d»fl?-\[A + ,A_} 2 + 

+ (d 2 a A_) A + d- 1 A + d~ 1 A + + (d 2 a A + ) A_dZ 1 A_dZ 1 A_ (259) 

and fl is determined by eq. (243). 

The term s 2 leads to the correctly normalized propagator 



(0 |T {A\{x))t{At{x'))i,\ 0) = 9(y - y' - rj) S 2 ^ - x[) ± In (* - x')\ (260) 

of the reggeon fields {A y ± )l (b, c are colour indices and A^,' = 8{8y — N^S^Sy is the 
projector to the adjoint representation of the colour group). Note, that we introduced y 
as an additional argument enumerating the Fourier components of the fields and included 
the factor 9(y — y' — if) in the Green function to remind that A± describe the interaction 
of the clusters with the relative rapidity bigger than rj. One can modify the term s 2 to 
guarantee this property: 

/+oo f) 
dyd^A y °d^A^. (261) 
-oo ay 

Here we took into account the interactions for all rapidities contrary to the case of the 
initial action defined in the fixed rapidity interval (y — rj, y + rj). 

The term S3 corresponds to the triple reggeon coupling. This reggeon vertex was 
introduced earlier by A. White [53]. Because the signature of the gluon is negative, s 3 
does not give any contribution to the elastic scattering amplitude due to the Gribov 
signature conservation rule but this vertex can be essential for some inclusive processes. 

The term S4 describes the four-linear interaction of the reggeons. The first two con- 
tributions correspond to their elastic scattering in the t -channel and last ones give the 
gluon transitions 1 <^ 3 which are not suppressed by the Gribov rule. The quantity /* 
is conserved: d a f^ = and is directly related with the effective RRP vertex C^gi,^)- 
Because for the big squared mass k ~d 2 of the intermediate state and fixed transverse 
momenta we have 

fldlfl^\[A + ,A„}\ 

there is a significant cancellation between two first terms in S4 leading to a better conver- 
gency of the corresponding integral: / where In k < rj. For the colorless state in the 
t-channel this integral is convergent in the ultraviolet region and by taking the residue 
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at k = one can obtain the result, corresponding to the real contribution to the BFKL 
kernel. 

The term s 5 also can be expressed in terms of the above calculated quantities f 1 and 
f 2 . The higher order term sq describes in particular the reggeon transition 2 — > 4 related 
with the triple pomeron vertex [54]. 

In a general case to cancel the infrared divergencies one should take into account 
apart from the classical expressions for the corresponding transition vertices also the 
contributions from quantum fluctuations near classical solutions. For example in LLA 
they are responsible for the gluon reggeization. 

Let us write down the field v as a sum of its classical component V and the small 
variation e describing its fluctuations near the classical solution: 

v = v + e (262) 

and expand the action in e: 

AS = S ef) - Sfff = - J d*xtr { [r>„,e„] 2 - [D„,e„][D„,e,J + gG w [e„, £ „] 

^-ib^-h^-^ + l^ik )(e+ ik mv+)A - + 0(£3) } ■ (263) 

The reggeon action in the semiclassical approximation can be obtained if one would calcu- 
late the functional integral over the quantum fluctuations e (with taking into account the 
Faddeev-Popov ghosts). Here we use the perturbative solution of the classical equations 
to write down AS only up to quadratic terms in e and bilinear in A±: 

AS = - J d 4 xtr {{d^f - {d^ u ){d u t^) + g{2(d /J ,e u )[v lx , e„] - 2(d 1/ e lx )[v^e l/ ] 
+2{d u v,)[e u ,e,} - (^ A+)e _J_ e _ _ (dl a A_)e + j-e + } 
+9 2 {[A + ,e„][A_,e u ] - \[A + , e + ][A_, e_] - \[A + ,e_] 2 - \[A^e + f 

+ (d{ a A_)(e + j-e + j-A + + e+ j-A + j-e + + A + j-e + j-e + )}^ , (264) 

where one should substitute v by the classical solution. The terms bilinear simultaneously 
in A± and in e can be used to find the next-to-leading corrections to the BFKL pomeron. 
In comparison with the dispersion approach [49] the method based on the effective action 
could provide a better infrared convergency of intermediate expressions. 

We consider below only the contributions which are linear in A± and bilinear in e CT 
with the singularities c^ 1 . Expanding the integrand exp (—iS e ff) up to the order con- 
taining linearly both A + and A_ and substituting the products of e + and e_ by the free 
propagators in the Feynman gauge 

f d 4 k 1 

< e_(x)e + (0) > = -2i J j—^ k2 + iQ expiikx) , (265) 
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we obtain effectively the box diagram because apart from the gluon propagators (k 2 + 
iO)^ 1 and ((q — k) 2 + iO)^ 1 there are also the propagators kz 1 and kz 1 appearing from 
the corresponding singular factors dz 1 and dz 1 in the vertices. The contribution of the 
Faddeev-Popov ghosts is small. As usual, one can obtain the logarithmic term ~ In s as a 
result of integration over k + and restricted by the constraint k + k_ ~ k\. We interpret 
it as the integral over the rapidity y of the t-channel gluons with momenta k and q — k and 
write down the corresponding one-loop contribution to the effective action in the form 

S] ff = - J d 4 x ± J d 2 x' ± J dytr(d ±a A y + (x.,x ± )) (d ±a A y S v (x + ,x' ± )) f3(x ± -x' ± ) (266) 

where /3(x±) is obtained by the Fourier transformation 

/d 2 q 
j^y ^(-^ 2 ) exp (i xl It) (267) 

from the gluon Regge trajectory: 

Note, that in LLA u(t) does not depend on the intermediate parameter 77 because it does 
not contain the ultraviolet divergencies in the rapidity. 

With taking into account apart from s 2 also the one loop correction the renormal- 
ized correlation function contains the Regge factor exp(u;(g 2 ) (y — y')) in the momentum 
representation: 

J d 2 x ± exp(2x ± g) < A y + (x)A y l(0) > ren ~ 9(y - y' - 77) exp(^(g 2 ) (y - y')) . (269) 

Analogously starting from the term ~(d^f u ) 2 in the contribution S4 (259) one can reduce 
the integrals over the momenta k± of the real intermediate gluon to the integral over its 
rapidity y and derive the corresponding term of the BFKL kernel. 

4.4 Effective vertices 

Because the classical extremum of the effective action is situated at non-zero fields v± = 
A± + ... , it is natural to parametrize v± as follows [51]: 

v ± = V ± + U(V T )A ± U~\V T ) , U(V ± ) = 1+ g i y± , (270) 

where V± is transformed under the gauge group similar to v± and V± = for small g. For 
this parametrization the expansion of S e ff in the series over A± is gauge invariant but 
the coefficients are rather complicated and some contributions to scattering amplitudes 
contain the simultaneous singularities in overlapping channels (cancelling after the use of 
equations of motion). It is the reason, why we shall use the simpler parametrization of 
v±: 

v± = V± + A ± . (271) 
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In this case one has also V± = at g = in the Landau gauge d a v a = and the 
homogeneous polynomials U of fields A± appearing in the expansion of S e ff\ 



/CO 
d A xtr £L\ 
i=0 



(272) 



are compatible with the Steinman relations. The terms L l do not have simple gauge 
properties but the corresponding scattering amplitudes are invariant under the gauge 
transformation after using equations of motion. 

In the above expansion L° describes the pure gluonic interaction: 



L° = ^(V) = - d v Vrf + 2g [V„ V„] 8,V V + 9 - [V M , V v ] 



(273) 



The next term L 1 contains the reggeon fields A± linearly and in particular it governs 
the quasi-elastic processes: 

L 1 = j+ A_ + j_ A + , 



(274) 



where the currents j± are given below: 



J± = -[D a ,G a± ]--dl a d ± ^d ± 
9 L> ± 



= 9 {[v a , d±V a ] - [d a V a , V±] - 2 [v a , d a v ± ] - dl v ± dz l v±) 
+g 2 ([v a , K, v ± ]] + d\ a v±d£v±d£v^ 

The quadratic term of the expansion 



(275) 



Dt~,A, 



D±,A. 



- \ ([D + , A.] - [D., A + ]f + | G + _ [A.,A+] 



+ ( d+— A + — d + - A + dl a A_ + ( d-—A_—d- - A_ ] tf„A 



D+ D 4 



D_ D. 



= d^A + d^A_ + g { V a \ A + \ ft A. + - L4_, A + ] d + V- - [V + , d~ 1 A + 



+ 



V a \ A_] d^A + + l - [A + , A_] d_V + - [v_, dl'A. 



d 2 ±a A, 



1 



1 



+g 2 1 [V a \ A + \ [V a \ A_\ — - ([VI, A + \ - [V + , A.])* + - [V + , V.\ [A., A + \ 

+ ({V + d- l fA + + A + (d- l V + r + V + d- l A + d- + l V + ) d\ a A_ 

+ ({V_dZ l ) 2 A_ + A_(9lV_) 2 + V_dZ 1 A_dZ 1 vJ) d\ a A + ) + ... (276) 

describes the gluon production due to the fusion of two reggeized gluons at the central 
rapidity region. 



46 



The following terms non-linear in fields A± are responsible for more complicated pro- 
cesses 

L 3 = 9 -([D + ,A_]-[D_,A + ])[A_,A + }- 

-I ( d +w + A +h A+ ^ d+dl ° A - + d -w^ A -wJ- d2 ^) 

= -| {A+d^A+^A. + A_dZ 1 A_d 2 ±a A + ) + 0(g 2 ) . (277) 
* = 4 [A + , A_f + £ (a, 1-A + 1-A + 1-A + ± d+dL A_ 

+d -w A -w A -w A -Tr d - d ^ A +) ■ ( 278 ) 

For i>4we have 

v = (-iy- y 




-1)! 

(279) 

The effective vertices obtained according to the Feynman rules from the above action 
are sums of the usual Yang-Mills couplings and some nonlocal induced terms. Below 
we construct the gluon production amplitudes in the quasi- mult i-Regge kinematics us- 
ing these effective vertices. One can find in the framework of this approach also the 
perturbative expansion of the reggeon action S regg defined as follows 

exp(-iS regg (A±)) = J DV eM-iSeff) (280) 

and depending on the reggeon fields A±. 

The subsequent functional integration over A± corresponds to the solution of the 
reggeon field theory defined in the two-dimensional impact parameter subspace with the 
rapidity playing the role of time. This theory is obtained after the integration of the multi- 
reggeon couplings over the reggeon longitudinal momenta within the fixed rapidity interval 
Ho — r], yo + f], which is equivalent to calculating certain limits of these couplings in the 
longitudinal subspace. It is important, that in the above approach the t-channel dynamics 
of the reggeon interactions turns out to be in the agreement with the s-channel unitarity 
of the S'-matrix in the initial Yang-Mills model. In the Hamiltonian formulation of this 
reggeon calculus the wave function will contain the components with an arbitrary number 
of reggeized gluons. Nevertheless, one can hope that at least some of the remarkable 
properties of the BFKL equation [46-48] will remain in the general case of the non- 
conserving number of reggeized gluons. 

Note, that to build the effective action for the multi-Regge kinematics one should take 
into account only two first terms of the perturbative expansion of L: 

L mR = \ (c^K - d u V,) 2 + d^A + d^A„ + gb 3 (A + , A_, V) , 
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h = ~\ A+d-'A+d\ a A_ - l - A_dZ 1 A_dl a A + + F + _ [A_, A + ] 

- fc^F+J) [d a A_,d a A + ] + (dz'F^) [A_,d a A + ] + (d^F^) [A + ,d a A_] 

-A+ [F-^F-v] - A_ [F^d^F+r] , (281) 
where we introduced the abelian strength tensor: 



F, u = d,V u - d u V, (282) 

and omitted in the last line some terms containing the factors d a V a and d 2 V^ vanishing for 
the real gluons. The Feynman vertices of this theory coincide with the effective reggeon- 
particle vertices of the leading logarithmic approximation. 

The obtained multi-Regge action is invariant under the abelian gauge transformations: 

V a ^V a + d a <j). (283) 

Its disadvantage is that the contributions of loop diagrams can contain the terms in- 
compatible with the Steinman relations forbidding the simultaneous singularities in the 
overlapping channels. Such unpleasant properties are absent for the full effective action 
(210) invariant under nonabelian gauge transformations. 

4.5 Gluon and quark pair production in the quasi-multi-Regge 
kinematics 

As an example of the general approach we consider the quasi-elastic process in which 
several gluons are produced with a fixed invariant mass in the fragmentation region of the 
initial gluon A provided that the momentum of the other particle B is almost conserved: 
p B i ~ p B [49, 51]. It is convenient to denote the colour indices of the produced gluons 
by di,o>2, ■■■O'n (a>i — 1,2, --N^ — 1) leaving the index a for the particle A. Further, the 
momenta of the produced gluons and of the particle A are denoted by k\, k 2 , ...k n and — k 

n 

correspondingly. The momentum transfer q = — h is fixed and t = q 2 . Omitting the 

i=0 

polarization vectors e Vi {k-i) for the gluons we can write the production amplitude related 
with the one-gluon exchange in the t-channel in the factorized form: 

Avom...v n . _ _±vovi-Vn+ 1 n rf-T Ax > (9RA) 

n -a a 1 ...a n B'B — <? 'a ai . . .a n c ^ 9 PB 1 B' B°^ b > ,A s • 

Here the form-factor <fi depends on the invariants constructed from the momenta ko, ...k n . 
In the simplest case of the elastic scattering the tensor equals 

4%£ = 9*Z£ (285) 
where after extracting the colour group generators 

Coalc = Kao ^^(kl: ~k ) (286) 



48 



the quantity T U1U0+ can be presented as the sum of two terms: 

r VlVo+ (ki, -k ) = Y lU0+ (h, -ho) + A™ + {k u -ho). (287) 
Here ^ U1U0+ is the light-cone component of the Yang-Mills vertex: 

Y lU0+ {h, -h) = S U1U0 (kt + A; + ) + {n + Y° {2k + k^ + (n + ) l/1 (-2^ - ko) u ° (288) 
and A UlUo+ is the induced vertex: 



A^° + ( M , -ko) = -t{n + y^{n + Y» , k+ + k+ = 0. (289) 

k{ 



For the case of the production of one extra gluon the amplitude was calculated 
several years ago [49] and in accordance with the above effective action (272) it can be 
written in the gauge invariant form 

Avoviv2+ _ 2 rp 1/0^1^2+ _ rpa rpc l Vll, ° u {ki-, —k ) T U2CT+ (k 2 , k 2 + q) 

Vaoaia2C iJ I 1 aoaia2C 1 a\a a 1 0,2,0, {^Zq + &i) 2 

Y 2U ° a {k2, -k ) T uicr+ (ki, k x + q) 



C 



- 1 a 2 a 1 1 o o 77, 7~i7Y2 / ' l 2yU J 



'■2 j 

r 2Uia (h, -fc!)r^ + ( fco, A; + g) 
(*i + A; 2 ) 2 

The last three terms in the brackets correspond to the Feynman diagram contributions 
constructed from the gluon propagator combining the usual Yang-Mills vertex 7 and the 
effective reggeon-gluon-gluon vertex T. The first term is the sum of two terms 

aoa\a2C I aoa\a2C ' aoaia2C i V / 

where 7 is the light-cone projection of the usual quadri-linear Yang-Mills vertex 

laoaia2C aiao H2iv u u u J 

+ K^TU^ 1 ^ - + T a a2ai T C a oa {^ Vl+ ~ 8^8^) (292) 

and A is the induced vertex appearing in the effective action (see (275)): 

KIZZUK, kf, = -t (n + r(n + r(n + r{%^ + 7 %^} . (293) 

ft l ft 2 ft ft 2 

Note that due to the Jacobi identity 

rpa rpc rpa rpc rpa rpc ( OQ/1 \ 

«2«o o\a a^a\ a§a a\OQ a^a v y / 

and the momentum conservation law 

kt + kt + kt = (295) 
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valid in the quasi-elastic kinematics the tensor A has the Bose symmetry with respect 
to the simultaneous transmutation of momenta, colour and Lorentz indices of the gluons 
0, 1, 2. The above expression for ^a a\a 2 2 t * s significantly simplified if we use the light-cone 
gauge e + (ki) = for the gluon polarization vectors because in this gauge all induced 
terms disappear. 

In a general case n > 2 for the gauge invariance of <p one should take into account apart 
from the usual Yang-Mills vertices 7 also an arbitrary number of the induced vertices A 
appearing in j + (275) and satisfying the reccurence relation: 

^oooi...o r c \ K ' K l ' — K r ) 

(n + ) Ur r ~ 1 

Ta r a i ^aoa 1 Lal^ 1 aa i+1 ...ar-ic(^0 > ■■■^t-li + > ^i+l> -"^r^l) ■ (296) 



k+ 



These induced vertices are invariant under arbitrary transmutations of indices i : 

Aa l ° a j L 1 ...ai r c (kf Q , kf v ■■■kf r ) = A^^""^ ( k£ , kf, ...kr) (297) 



due to the Jacobi identity for the colour group generators T and the energy-momentum 
conservation 

XX = 0. (298) 

i=0 

One can calculate easily also amplitudes of the quark production in the quasi-elastic 
kinematics because the reggeized and usual gluons interact with quarks in a similar way. 

Let us consider now the multi-gluon production in the central rapidity region [49, 51]. 
The following kinematics of the final state particles is essential: the gluons A' and B' move 
almost along the momenta of the initial gluons A and B and there is a group of produced 
gluons with a fixed invariant mass in the central rapidity region: y = |ln(A; + /A;~) ~ 0. 
The momentum transfers qi = Pa — Pa> and q 2 = Pb' — Pb in this regime have the 
decomposition 

Qi = Qi± + Ppa , <?2 = <?2_l - aPB (299) 

n 

where (3 and a are the Sudakov parameters of the total momentum k — J2 &i of the 

i=i 

produced gluons: 

k = k± + f3p A + aps , k = k 2 = sa/3 + (qi - q 2 )\ (300) 

and y/R is their invariant mass which is assumed to be fixed at high energies: k<s. In 
this kinematical region the production amplitude has the factorized form: 

1 1 

/i d 1 d 2 ....d n A'AB'B — ~ 9Va 1 A'A°^a^ a i T~ t f J d 1 d2...d n c 2 c 1 T~ 9 Pb 1 B' B °A B A B , • 

For the simplest case of one gluon emission we have 

Vd\i^=gT7Xcv ( 302 ) 



50 



where T is the sum of two terms 

^dlt 2 g = IdlXci + ^d\t 2 ci- (303) 

The first term is the contribution from the tri-linear Yang-Mills vertex 

7g+- = T^ + ~ , 7 = 2(q 2 + qi ) - 2kfn~ + 2k^n + . (304) 
The second term is the induced one 



71 71 

K\tg = TtA V1+ ~ > A = "2 *i F + 2 1 2 tt • (305) 



Due to the relation 

7 + A = -2 C, (306) 

we obtain the known result for the multi-Regge kinematics. 

For the more complicated case of the two gluon production in the central rapidity 
region the amplitude was calculated also and the quantity ip is given below [49, 51]: 



Yd 1 d 2 c 2 g y I 1 d 1 d 2 c 2 c 1 ^ _|_ fc^2 



(qi - h) 2 (qi - k 



\2 



2) 

The second term in the brackets describes the production of a pair of gluons through 
the decay of the virtual gluon in the direct channel. This contribution is a product of the 
effective vertex T, the usual YM vertex 7 and the gluon propagator. Analogously, the 
third and fourth contributions are products of two effective vertices T having the light 
cone components ± and of the gluon propagator in the crossing channels. 

The first term in the brackets can be presented as the sum of two terms 

-pviu 2 A — _ -\>\v 2 -\ — 1 \v\v 2 -\ — fQn8"\ 
1 d x d 2 c 2 g — ld 1 d 2 c 2 g L *d 1 d 2 c 2 gi [JUS) 

where the contribution 7 is the light cone component of the quadri-linear Yang-Mills 

vertex 

+ Tl cl T£ d (8™8- + - + T d d2d XU5 V2 -^ + ~ (309) 

and A is the new induced vertex appearing in the effective action (see (276)): 

rpd rpC 2 rpd rpC 2 

\v\v 2 ^ — _ 9/ L + ^L+^/fclW 1 d 2 di gd -| 

A did 2 c 2 ci - -*t 2 \n ) [n ) 1 7+ 7+ + 7— r+ — , +] r+i ~ 

rpd rpg rpd rpg 

- 2ti (n-p(n-y 2 { dl ? d J d + dld2 c l d — }. (310) 

k 2 fei {~ki — k 2 ) /e 1 

One can verify that A is symmetric under simultaneous transmutations of n~,ti,di and 
n + ,t 2 ,d 2 . 
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Note, that in the general case of multi-gluon production one should take into account 
the induced vertices with an arbitrary number of external legs which are expressed in 
terms of the light cone projections of the vertices introduced above for the quasi-elastic 

case 

^dl'.'.'^Xci = ^tid\''. U d n c 2 (^0^1---^n) + ^SV-.^ci (^0 ^1 -■^n ) > (3H) 

where k£, k$ are determined by the momentum conservation 

n n 

E*+ = E*r = o. (312) 

i=0 i=0 

The tensor ipdld^cx appearing in the two gluon production amplitude A 2 -,a can be 
written in the other form: 

C*£« = 29 2 [TdldTZd^ihM) +TZ d T? ld A™\k 2 M)] , (313) 

where 

r ai/3 ~(fci, h - qi ) T a2 ^(k 2 , k 2 + q 2 ) r 2ail3 (h, ~h) T+-P(q 2 , Ql ) 



2 (qi - 2 (h + k 2 f 



+ n +ai n~ Q2 - g aia2 - - n +a2 n - ai + t^—— — + t 2 -—+ — . (314) 

2 /c^ [k-y ~~t~ k 2 ) k 2 [k-y ~~t~ k 2 ) 



One can verify the following gauge property of A" 1 " 2 : 

L«u«ia 2 _ l u*2 ( fe i k 2 fcfr + l3 (q 2 ,qi) \ , . 

^ " 2^ 2 il^T^F (fci + fc 2 ) 2 J ' (3 5) 

which implies, that the production amplitude multiplied by the gluon polarization 

vectors e{kj) is gauge invariant. To calculate the matrix element of A ai ° 2 it is convenient 
to use the different light-cone gauges for two produced gluons: 

e~(k 1 ) = 0, e + (k 2 ) = 0. (316) 

The corresponding polarization vectors can be parametrized by the two-dimensional vec- 
tors eij_, e 2 ± as follows 

a \ (foei±) _ n v (^2e 2 ±) + _v 
e(«i) = ei ± — — n , e[k 2 ) = e 2± T^ n (317) 

K-y K 2 

and the matrix element of the tensor A aia2 can be expressed in terms of a new tensor 
c" 1 " 2 with pure transverse components according to the definition: 

^^iX.W^^e^^j^Wc^. (318) 

Note, that the corresponding tensor Z a2ai for A a2ai (k 2 ,k 1 ) is obtained from c by the 
transmutation of gluons 1 2 with the simultaneous multiplication by the matrices 

k ai k ai ' 
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interchanging two light-cone gauges (316). 

Using also the reality condition k + k~ = \ki\ 2 for the produced gluons i = 1,2 to 
exclude the light-cone momenta k~ one can present the tensor c" 1 " 2 only in terms of 

transverse momenta ki , ~qt and the relative Sudakov parameter x = k +^_ k + [55]: 



^ lQ2 _ 2 , t x , ai t2 2 K 2 + ^ 2 QT - (i - x)-\qT - t 2 )k? 



J_ , QinQ2 , A°i£>+6™qI(K + x<8)-(l-x)- 1 £ 1 (k?-xA<») 



K 



QTQV - |(1 - fci )5 ai " 2 2 ^2 ^2 2 + 5 aia2 (l — A 
+ ^ , (319) 

where 5 aia2 = — g" 1 " 2 is the Kroniker tensor (5 11 = 5 22 = 1) and 

A = 5i - q 2 = h + k 2 , Qi = qi - h , k = — , 

x(l — x) 

Z = -X(l -x)(Z 2 + K),t= _ (^-^) 2 +*(l-^ 2 . (320) 

X 

Analogously to the gluon case, the amplitude of the production of the quark and anti- 
quark with their momenta k\ and k 2 in the central rapidity region also can be written as 
a sum of two terms being the matrices in the spin and colour spaces [55]: 

^ C2C1 = -g 2 (t c H C2 b(k u k 2 ) - t c n c ^b T (k 2 , h)) . (321) 

Here t c are the colour group generators in the fundamental representation and the expres- 
sions for b(ki,k 2 ) and b T (ki,k 2 ) are constructed according to the Feynman rules including 
the effective vertices: 

uu u\ + 7^r + - /3 (g 2 , qi) fQOO , 

b{h • h) = 1 w^r- ' - ( kl + k 2 f- • (322) 

urn / \ + qi-fo _ 7 /3 r + - /3 (g 2 ,gi) , QOQ ^ 

* *.^) = •» - • (323) 

One can calculate the matrix elements of these expressions between the spinors u^(ki) 
and v^(k 2 ) describing the quark and anti-quark with helicities A = ±|: 

^(hWhukdv^fa) = ^c* + _(k-i)d- + (k 2 ) b^ + ~\h, k 2 ) , 

u^ih^ik^kjv^ih) = ^c*_ + {k 1 )d + -{k 2 )b ( '- + \k 2 ,k 1 ), 

u^\ki)b(ki, k 2 )v i+ \k 2 ) = ^(kjd—ih) b^ + \h, k 2 ) , 
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u { -\h)b T (k 2 , h)v( + \k 2 ) = ^c*_Uh)d ++ (k 2 ) b^ + -\k 2 , h) . (324) 
The introduced complex wave functions c and b satisfy the relations 



k\ k\ . . k 2 _ /c : 



c+-(fci) = -^c_+ = -f^c_+ , d + _(/c 2 ) = -tt^-+ = - jr d -+ , 



h* 1c h* h 

A/^ rv 2 rv> 



2 



c—(ki) = ~jpc ++ = -^-c++ , d—(k 2 ) = ~^d ++ = ~j^d ++ , (325) 
and are normalized as follows 

| c&ih) | 2 = 2k* , | d ±l (k 2 ) \ 2 = 2k% . (326) 
in accordance with the usual normalization condition for the spinors: 

uik^uih) = 2k? , v(k 2 )-f a v(k 2 ) = 2k% . (327) 
The amplitudes b can be written in the form 



,(+-)/, ,n /,,(-+) /i 7 \ \* , /r+ "(g 2 ,gi) 

6 ) (^^) = ( 6( )( ^ fc2)) = - 4 (^^~ (fcl + A;2)2 



where the quark current j is given below 



(328) 



. * k 2 k\ _ k\ + k 2 /oon\ 
j = n + n — ra rv — . (329) 

rb 2 ?h\ rh\ rv 2 

Due to the relations 

n a k a = k, n* a k a = k* , n ±a k a = fc* 

this current is conserved: 

(*f + *£)/ = 0. (330) 

5 Next-to-leading corrections to the BFKL pomeron 

5.1 Corrections from the gluon and quark pair production 

The imaginary part of the elastic scattering amplitude calculated with the use of the 
s-channel unitarity condition through the squared production amplitude A 2 ^ (301, 313) 
contains the infrared divergences at small kj ± and k. To avoid such divergencies the 
dimensional regularization 

^-d 1A (331) 



(2tt) 4 ^ (2n) D 

is used in the gauge theories (where \i is the normalization point). It is important, that 
in the .D-dimensional space the gluon has D — 2 degrees of freedom. 
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The total cross-section of the gluon production is proportional to the integral from 
the squared amplitude V'di^ctcT summed over all indices. It is expressed in terms of the 
bilinear combinations of the tensor c" 1 " 2 with transverse components [55]: 



R = (c aia2 (ifei, k 2 )f + {c aia2 {k 2 , k x )f + c aia2 {k u k 2 ) c a '^{k 2 , h) Q aia 'i Q a2a *, (332) 
where the matrix VL ai0ii ' interchanges the left and right gauges: 

no.",' _ „ a i a i> oVrii 

" 9± ~ p — • 

K i± 

In expression (332) we took into account that the colour factor for the interference term is 
two times smaller than one for the direct contributions. The generalized BFKL equation 
describing the pomeron as a compound state of two reggeized gluons is valid also in the 
next-to-leading approximation. The equation for the total cross-section of the reggeon 
scattering can be written in the integral form as follows 



/d q^ f 
/i 4 ~ D J (fql K 5 (qi, q+ ) . (333) 



Here the integration region for the longitudinal momentum q 2 is restricted from above 
by the value proportional to q x : 

qt < 8q+. (334) 

The intermediate infinitesimal parameter 8 > is introduced instead of the above cut-off r\ 
to arrange the particles in the groups with strongly different rapidities. The integral kernel 
KsiqXiQ?) takes into account the interactions among the particles with approximately 
equal rapidities. The kernel K$ can be calculated in the perturbation theory: 

oo / 2 \ r 

We remind, that the real contribution to the kernel in the leading logarithmic approxi- 
mation is proportional to the square of the effective vertex: 

K%M,-fi) = -7^2 (r^te,*?!)) 2 = N c 4 (336) 
4gi g 2 v 1 (<?i - Q2) 2 

The corresponding virtual contribution is determined by the gluon Regge trajectory 
The next-to-leading term in Kg related with the two gluon production is given below 

[55] 

K W f dK fh S(kW) ^SE 16 N < r dx f dD ^ R (337) 
K giuons - J an ^ D _ 4 8{k 1 ) 8{k 2 ) ^2-^2 - 2 ^ 2 ^ 2 J s x(l - x) J {667) 

The limits in the integral over x correspond to a restriction from above for the invariant 
mass yfk of the produced gluons. In the solution of the generalized BFKL equation the 
dependence from 8 should disappear. 
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For the physical value D = 4 of the space-time dimensions one can express R as the 
following sum: 

R = R(+-) + R(++) (338) 

where R(-\ — ) and R(++) are the contributions from the production of two gluons with 
the same and opposite helicities correspondingly. 

R(+~) = l(\ c + -(h,k 2 ) | 2 + | c + -(k 2 ,h) | 2 +Rec + -(k 1 ,k 2 )c- + (k 2 ,k 1 )^-^j , 

R(++) = \(\ c ++ (h,k 2 ) | 2 + | c ++ (k 2 ,h) | 2 +Rec ++ (k 1 ,k 2 )c ++ (k 2 ,k 1 )^ 1 ^j . (339) 

Here we have used the following relation between the polarization vectors (317) in the 
right and left gauges [42]: 

ei(fc) = (340) 

to express amplitudes (319) in terms of two complex functions c + ~{ki, k 2 ) and c ++ (/ci, k 2 ) 
given below: 



c + -(k u k 2 ) = c 11 + ic 21 - ic 12 + c 22 = c-+(h, k 2 ) 



-x 



(ki -xA)kf 

c ++ (/ci, k 2 ) = c 11 + ic 21 + ic 12 - c 22 = 7r^(ki, k 2 ) 
x(Qi) 2 + xqt 2 (k 2 ) 2 



((£ - xqlf + x(l - x)ql 2 ) ^ 2 ^ _ x ~£y + x (i _ x ^ 2 ^j 

_ x(l - x)q l k 2 g* 2 _ xq\k x q 2 xq 2 Qi 
~ A*(h-xA)kt ^\ki- x A*) A* kf 

These expressions were obtained independently also in Ref. [56]. 

In the Regge regime of small 1 — x and fixed k iy q^ amplitudes c aia ' 2 are simplified as 
follows 

c + (h, k 2 )-^-—, c ++ (h, k 2 ) — > — — — — (342) 

and are proportional to the product of the effective vertices r +_/3 in the light-cone gauge 
[42]. For x — > the amplitudes c + ~(ki, k 2 ) and c ++ (k\, k 2 ) vanish, but for simultaneously 
small ki or k-y — xA one obtains a nonzero integral contribution because in this region 
there are poles: 

For large k\ and fixed qi , x we obtain 

c + (/ci,ft 2 ) — ► -x— 2 , c ++ (fci,fc 2 ) — ►-a;(l-a;) — T - a;' 3 — — - (344) 
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and therefore the integrals for the cross section of the gluon production do not contain 
any ultraviolet divergency. 

As it is seen from above formulas, c aia2 has infrared poles at small k\\ 



•(*.,*,)- ff.^.M-ff (345) 



and at small k\ — x A : 



S^-V-x? T^P^ -J ■ (346) 



K {2) 

quarks 



(347) 

Jo x(l -x) J at at 



A*(h-xA) y ' (A*) 2 (k 1 - xA) " A*(kl-xA* 

Let us return now to the quark- ant iquark production. The total cross-section of this 
process in accordance with the normalization condition for spinors is proportional to the 
integral from the expression: 

- 1 dx^' D f jD _ 2T > 8L 

where we put 5 = because the integral in x is convergent at x = and x — 1. 
The expression for L with corresponding colour factors for D = 4 is given below: 

/V 2 — 1/1— r r \ 1 

L = [—— | c(h, k 2 ) | 2 +— | c(k 2 , h) | 2 j + —Re c(k u k 2 ) c(k 2 , h) (348) 

where two equal contributions from two different helicity states of the quark- ant i- quark 
pair were taken into account. The amplitude c is determined by the equation (see (328)): 

b (+ -\k 1 ,k 2 ) = ^c(k u k 2 ). (349) 



and is given below: 



1 (1 - x)q 1 k{ - xq^h + xq?- 

-c(k 1: k 2 ) - 



x ' (fci - xqt) 2 + x(l - x)qt' 



qt 2 (1 -x)Ak{ -xA*h + xA~ (l-x)q l q* 2 xq^q 2 



A*' 2 (kl-xty + x(l~x0 2 ' A*(h-xA) A(kl - xA*) ' (350) 

Note, that in the limit of large N c the quark contribution is smaller than gluon one 
and the interference term is suppressed by the factor 1/iV 2 in comparison with the direct 
contribution. 



5.2 Infrared and collinear divergencies 

The gluon and quark production cross-sections contain infrared divergencies which should 
be cancelled with the virtual corrections to the multi-Regge processes. For example from 
the products of amplitudes c + ~ we can calculate R(-\ — ) for D = 4 and obtain the following 
expression: 

R ( + _) = ^ 2 ( 1 + x " + t 1 - x ? \ ( 351 ) 

4 V^ 2 (^-^) 2 kl 2 (M-xA*y (kl-A*nkl-xA*y) 
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from which the presence of infrared and collinear singularities is obvious. 

As for R(++), we can write it as the sum of singular and regular terms [55] 

R(++) = R sing ( + + ) + R reg ( + + ) (352) 

where for R S i ng (++) the following expression is chosen: 

Rsin 9 (++) = R(+~) + \ (r(h, x) + r(A — ki,l — x)) . (353) 

and 

/, v giV , n x ^ 2x 2 {l -xjkTZ op x 3 (l-x) 2 q 2 q* 2 A 
d/ '-'- r) {X(1 ~ X) - 2) t\t- X ^) 2 + 2 ^ A* 2 ( kl -xA) 2 kl • 

Therefore for the total sum R S i ng = R{~\ — ) + R S mg(++) we have 

R sing = 2 R(+-) + 1 (r(A;i, x) + r(A - fci, 1 - x)) . (355) 

The term R S i ng contains all infrared and Regge singularities of R and is generalized to the 
case of arbitrary D. It decreases rapidly at large k\. 

One can obtain the following contribution from R S i ng after its dimensional regulariza- 
tion (and taking into account D — 4 non-physical polarizations for each gluon): 



+ 



^2 T( D^3 )sin(7r D-^ 



/i 2 



(z^4 +41 4-Y)' ^ 



2 / V" 2 

where it is implied, that the terms of the order of value of D — 4 should be omitted. 
The infrared divergencies at D — > 4 in the above formulas should be cancelled with the 
contribution from one-loop corrections to the Reggeon-Reggeon-particle vertex which will 
be considered below. 

It is important to investigate the region of small A: 

|A| < \ qi \ ~ \q 2 \ . (357) 

Here the essential integration region corresponds to the soft gluon transverse momenta: 

fci ~ k 2 ~ A < g (358) 

where g means gi or g 2 . 

The total contribution of the singular and regular terms in the soft region A — > can 
be calculated [55]: 



r 3 — — ^ — fx 
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D-2 

+ 



i) + 2 ^w^ ( 4 Xn \ + + ^ - l) - m ° - 3) ) 



(D-4) 2 2(D-4)(D- 1) 

(359) 

In this expression three first terms of the expansion in e = D — 4 including one of the order 
of e are correct. Such accuracy of calculations is needed because the integration over A in 
the generalized BFKL equation leads to the infrared divergency at A = for D — > 4. It 
means, that to obtain the final result after cancellations of the infrared divergencies one 
should know the real and virtual corrections to the BFKL kernel up to terms ~ e. 

Note, that it is possible to modify the definition of the singular and regular parts of 

c ++ (h,k 2 y. 

c ++ {h, k 2 ) = 4+ (h, k 2 ) + c++(h, k 2 ) , 
1 ~++ n / \ ql&hih - xA) - (2 - x)q x q* 2 Ak 1 + q^* 2 A 2 

~ C sing\ K l' K 1) - 

Jb 



A*' ((kl - x^) 2 + x(l - x)^ 



>2 ql<l2kl + (2 - x)^g - (1 - xf 



~&(kt-xA*) A * k *i A^h-xA)^ 

1 &h,k 2 ) = - Qf 



x reg {k\ — xql ) 2 + x{l — x)qi 2 

, qlhjh - xA) - (2 - x)q 1 A*k 1 + gi A*' q*h 



((kl - xt) 2 + x(l - x)t 2 ^j A * k ± 



in such way to include in the comparatively simple expression for c^ g (ki, k 2 ) all singular 
terms without the loss of its good behaviour at large k±. In the Regge limit x — > 1 and 
fixed ki we obtain 

and therefore in this region the contributions of the exact amplitude c ++ and approximate 
one cj"in S to the differential cross-section also coincide. However, as a result of the singu- 
larity at k\ = xA after integration over k± these contributions to the total cross-section 
in the Regge limit x — > 1 turn out to be different. 

We return now to the quark production amplitude c(ki, k 2 ) and present it in the form 
analogous to the gluon case: 

c(h, k 2 ) = Csingih, k 2 ) + c reg (k!, k 2 ) , (362) 

where the singular and regular terms are chosen as follows 



— >2 

1- ( , , s _ (1 - x)q 1 q 2 xq{q 2 (1 - x)Aq 1 q 2 f k* 1 - xA*q* l q 2 k 1 + xA q^q 2 
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1 (1 - z)<7i£;* - xq{ki + xg^ 2 (1 - x)q 1 k\ - xg*A;i + xg*A 

x ^k 1 -xq i y + x{l-x)q 1 ) h kl - x A) 2 + x(l - x) A J 

The term c refl does not lead to any divergency and gives a regular contribution to the 
BFKL kernel at the soft region A — > 0. 

We obtain the following contributions for the bilinear combinations of c S i ng taking into 
account the dimensional regularization of the singularity at k\ ->iA [55]: 



- 1 dx r d D 2 ki 2 f 1 dx f d D 2 k\ 

J a D ~H l c ^^i^2)l =y o _ x ) y „D- 47r Recsingih, k 2 )c sing {k 2 , fci) 



l /.TV \. . n r(#)r(£) . /6-£> 



^ 2 ^ 2 y ^"i^ 4 [d=Z* » ») J • ( 364 ) 

The integration over A in the generalized BFKL equation leads to the infrared diver- 
gency at A = for D — > 4: 

The quark correction to the RRP vertex expressed in terms of the bare coupling constant 
does not give such singular contribution at small A . This divergency is cancelled with 
the doubled quark contribution to the gluon Regge trajectory (see below). 



5.3 One-loop corrections to reggeon-particle-particle vertices 

To find next-to-leading corrections to the BFKL equation one should calculate also the 
analogous corrections to the gluon production amplitudes in the multi-Regge kinematics 
[41, 49]. These amplitudes in LLA are constructed from the RPP and RRP vertices 
combined by the reggeized gluon propagators expressed in terms of the Regge trajectory. 
In the next-to-leading approximation this structure of the production amplitude remains 
and therefore one should calculate one-loop contributions to the reggeon-particle vertices 
and two-loop corrections to the gluon Regge trajectory. 

To begin with, we consider the elastic amplitude for the gluon-gluon scattering. In 
the Regge model it can be written in the form: 

T c B/B r XBlXB (t), t = q\ (366) 
where T\ A ,\ A {t) due to the parity conservation can be presented as the sum of two terms: 

r XAlXA (t) = r«(t) s Xa „ Xa + r( 2 )(t) 5 Xa „_ Xa (367) 

in the helicity representation or 



M = T% A T XAlXA {t)- t 



-) + (- 
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iW*) = r«(0 (-^) + r( 2 )(t) (^i - 2^) (368) 

in the tensor representation. Here 

" #i = E 4(PaO #(pa) , e a p« = 0. (369) 

A=± 

is the projector to the physical gluon states in the light-cone gauge. The vectors and 
q^i are orthogonal to the momenta k, Pb and k', Pb correspondingly . For finding the 
next-to-leading correction to it is enough to calculate in the one-loop approximation 
the scattering amplitude at high energies. The t-channel unitarity is convenient for this 
purpose. In the framework of this approach the t-channel imaginary part of the gluon- 
gluon scattering amplitude in the pure gluodynamics is obtained as a product of two 
Born amplitudes with summing over the colour and Lorentz indices and integrating over 
momenta —k, q + k of the intermediate gluons in the octet state. Each of these amplitudes 
is presented as a sum of two terms after projecting to the octet state. For example, 

A^'ci(k,q + k,p A )=g 2 (-n, A T d BIB ) (A& + A™)* , (370) 

where A^ is the factorized contribution containing the correct asymptotic behaviour of 
the total amplitude: 

4^ = 6£A*(k,k + q,p A ) (371) 

and A^ is the non-asymptotic term which is small in the Regge kinematics. The factor 
A 11 (k, k',p) is proportional to the effective RPP vertex: 

^ ( ^, p) = _ 2 ^^ + 2 (_i_ + |) J) ^ 
+ 2 ((F^p + l) ^ + 2 (wTW ~ whr) pV - k - k + q - (372) 

This tensor has the simple gauge properties on the mass-shell k 2 = k' 2 = 0: 

A 11 ' (k, k', p)k J = 0, A 11 ' (k, k', p) k' 1 ' = (373) 

and therefore the Faddeev-Popov ghosts are absent in the t-channel imaginary part of 
the one-loop amplitude. Further, the product of two non- asymptotic terms A^ does not 
give any contribution of the order of s to the elastic scattering amplitude. The integral 
contribution from the product of the terms A^ is given below in the D-dimensional 
space-time [41]: 



„4 

A(aa)aa'PI3' _U 4 ™«' r/3/3' a t T d T c 
71 AA'BB' — ~r9 °±±°±± ly c-l-A'A 1 B'l 



r(2-D/2) [T(D/2-l)] 2 

~ X 



(_t)2-r>/2 (4n) D / 2 T(D -2) 
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' ' - ' l + (D-3)( In Lil£ + 2 ^(3 - D/2) - A^{D/2 - 2) + 2^(1)' 



D - 1 D - 4 2 v y \ t 2 

(374) 

In an analogous way the contribution from the interference term between and 
AW equals [41] 



Ai na)a*>W g 4 ^ d 2 s T(2 - D/2) [T(D/2 - 1)] 



t ) 2 -D/2 (4vr)^/ 2 r( J D-2) 



■X 



1 2 

+ 



D-\ D -4 



+ 



t L> - 1 



Using the following regularized expression for the gluon Regge trajectory 



(375) 



(376) 



(4tt)^/ 2V r(£>-3) 

we obtain for the unrenormalized form-factors T^(t) and r^ 2 ^(t) of the reggeon coupling 
with gluons [41]: 



r (1 \t)=g + g" 



N c T(2 - D/2) [T(D/2 - l)f 



x 



(_t)2-D/ 2 (4ir) D / 2 T(D-2) 

7 1/4 



(D - 3) (^(3 - D/2) - 24>(D/2 - 2) + ^(1)) 



and 



r( 2 )(t) 



D - 1 

iV c r(3-D/2) [T(D/2- l)] 2 1 



(377) 



(378) 



[-t) 2 - D / 2 {4n) D / 2 F(D - 2) D - 1 

Therefore in the pure gluodynamics in the one-loop approximation the s -channel helicity 
of colliding particles is not conserved even in the limit D — > 4. The contribution to the 
reggeon-gluon-gluon vertex from the virtual quarks was calculated in Ref [57]. For rif 
massless quarks the result is 



AT^(t)=2n f g 3 (-t) 



D/2 _ 2 r(2 - D/2) [T{D/2)f 
(4tc) d / 2 T(D) 



At [t)~ 2n f g { t) (4tt)^ T(D) 



(379) 



For D — > 4 we obtain for A T (2) : 



Ar' 2 »(^r( 2 >W(-^). 

The total contribution of quarks and gluons to the vertex T^ 2 ^ responsible for the helicity 
non-conservation is zero for rif — N c — 3 . This fact is a consequence of the super- 
symmetry. Indeed, in the case of the super-symmetric Yang-Mills theory the helicity of 
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the gluon is conserved because it takes place for the gluino belonging to the same super- 
multiplet. On the other hand, this conservation in the one- loop approximation appears 
as a result of the cancellation between contributions of the virtual gluon and gluino. The 
gluino contribution can be obtained from quark one by its multiplication by the factor 
N c /rif because the gluino belongs to the adjoint representation and it is invariant under 
the charge conjugation. It leads to the vanishing of the helicity non-conserving transition 
at rif — N c — 3. 

The one- loop correction to the reggeon- quark- quark vertex was calculated in ref. [57]. 
In the massless case we have: 



r q (t)=g + g 3 (-f) 



D/2 _ 2 T{2-D/2) [Y{D/2-l)X 
{A<k) d / 2 T{D-2) 



n_ f 2 
2 



D 1/3 D 



D — 1 N r 



D-A 



+N C 



(£>-3) (v(3-y)-2V»(y-2)+^(l 



1/4 



D — 1 D-A 



(380) 



The infrared and ultraviolet divergencies can be extracted easily from expressions (317- 
380) and are in the agreement with the general considerations. 



5.4 Loop corrections to the reggeon-reggeon-particle vertex 

The amplitude M 2 ^ for the production of a gluon in the central rapidity region can be 
written in the Regge model as follows [41] 



M. 



2^3 



*n>A (f i) r T ti T T b>b Ta s ,a s (t 2 ) f(s, si, s 2 , qt, tf) , 



tl 



1 

h 



/-J 



+ 



W1—L02 



[s^ + i-sf 2 ] R(t u t 2 , k ± ) 



LU 2 -CUl 

b 2 



+ (-82) 



[s^ + i-sf 1 ] L(t!,t 2 , k 



±J 



(381) 



where j\,2 = 1 + ^1,2 (for iO{ = u>(ti)) are the gluon Regge trajectories with the negative 
signature at the crossing channels ti,t 2 . It is important, that this amplitude does not 
have the simultaneous singularities in the overlapping channels si = 2kpA and s 2 = 2kps- 
In LLA one obtains: 

R(t u t 2 ,t 2 ± ) + L^t^tl) -> 2gC^q 2 , qi )el(k) (382) 

where is the RRP vertex. One can calculate the discontinuities of M 2 ^ 3 in the s, si 
and s 2 channels for small g with the use of the unitarity condition in the one-loop ap- 
proximation in the D-dimensional space: 



R-L 



C^q 2 , qi )e*(k) AN c g 



D 



L0 1 - LV 2 



( 47r )D/2 r ( 3 2 



In 



1 



k 



D-A 



(383) 



Note, that this quantity is real in all physical channels. By adding the analogous contribu- 
tion of the order of g 3 to the sum R + L in such way to provide correct analytic properties 
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of the production amplitude in s, Si and s 2 channels we obtain for the one-loop correction 
to M 2 ^ 3 [41]: 



M l °Z 



M Born N - F(3 — —)- 

M 2^3 *c {4n)D/2 L{6 2 ) 2 



A-D 



In 



;- s ) 3 sV 



1 A„ 2 (si)(s 2 ) 



--lln 



—s 



+ ln 2 



SlS 2 , ,2 Sl(-S 2 ) , , 2 (-Sl)«2 



— s 



+ In 



+ ln 2 



+ AM 2 ^ 3 , (384) 



where AM 2 _> 3 is an polynomial in A;^ = Sis 2 /s. One can reproduce this analytic struc- 
ture of M 2 _+3 using the approach based on the unitarity conditions in t\ and t 2 channels 
and obtain the following result [41]: 



R + L = 2gC ll (q 2 ,q 1 )e;(k)f 1 + 



Pb_ _ Pa\ 

s 2 s l ) n 

h ~ l v 1 + Cl (47p^ r(3_ ¥V' /2 " C2 (4^ 

Ci = - In 2 k 1 - 

2 ± 



e;(k)Agf 2 , 



r(3-f); 



1 h + t 2 



3 (t! - t 2 ) ; 



+ d^a (y + 2 ln {tlh) ) + \ ln2(M2) 

6 * 2 V *i-*2 (ti-t 2 yj 2 



C2 = Iln^fll_ 



2 2 fc 



1 1 + f 9 



tit 



lf-2 



t 2 



(tl " t 2 ) 2 ) ti - t 2 



tl+ ' 2 (2tl + t 1+ t 2 )V 



(ti-t 2 ) 2 



(385) 

The contribution to the RRP vertex from the massless quark loop was also calculated 
[58]: 



N c 2/3 ~t\ h + t 2 1 ^ /*! + 1 2 oT >2 t x t 2 \ iV c 

^ AC1 = Z^4-— (^F + 3 ln ^ i^ + 2 ^(^Fj' 



>2 \ — >2 



=iin^ - 2 -t: 2 ^ +t \t t2 1+ ji±\ (2 t:+t 1+ t 2) 

3 t 2 ^ X (ti-t 2 ) 2 )t x -t 2 Q \ (t!-t 2 ) 21 x 2 V 

(386) 

and therefore for the super-symmetric Yang-Mills (where effectively rif — N c ) the total 
gluon and quark contribution to c 2 is significantly simplified. 



5.5 Two-loop corrections to the gluon Regge trajectory 

For finding the next-to-leading corrections to the gluon Regge trajectory one can use the 
unitarity relations in the t or s channels. Initially the t-channel unitarity approach was 
advocated [49]. In the framework of this approach the universality of the Regge trajectory 
obtained from the different high energy processes is obvious. To find the next-to-leading 
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correction to the gluon trajectory one should calculate the imaginary part of the elastic 
scattering amplitude in the two-loop approximation taking into account the terms of the 
order of ln 2 s and Ins. The contribution from the two-particle intermediate state in the 
t-channel can be expressed in terms of the product of two elastic amplitudes. One of 
these amplitudes is taken in the Born approximation and another one corresponds to the 
helicity conserving part of the above calculated one-loop elastic amplitude in the Regge 
kinematics. As for the contribution from the three-particle intermediate state in the t- 
channel, one should take the product of two amplitudes A 2 ^ for the quasi-multi-Regge 
kinematics for each of two inelastic amplitudes A 2 ^. Because in the leading logarithmic 
approximation the Regge expression for the scattering amplitude has only two-particle 
intermediate states, it is enough to calculate only the asymptotic behaviour of this three- 
particle contribution. 

One can use also the s-channel unitarity to find the asymptotic behaviour of the 
scattering amplitude in the two-loop approximation [59]. In this case the imaginary parts 
of the amplitude in the s and u channels should be calculated with taking into account 
the terms of the order of s In s and s. Subtracting from the result the known one- 
loop corrections to the reggeon-particle-particle vertices, one can calculate the two-loop 
corrections to the Regge trajectory. For the case of the gluodynamics it is given below 
[59]: 



N 2 c gH r d D - 2 qi r d D ~ 2 q 2 
4(27r) 2 ^- 2 J ~W~ 



— >2 — ¥2 o — >2 

, q , q 2 i ft 

J = —— 2 —2 In — : — -2 + — — — -J In — — -j 

(9-92) (q - qi) (?2 - qi) (q - qi - 92) ( 9 - 9i) 

>2 



+ 



q 



'cf-qi-qZ) 2 (^-92) 2 (g" 




(*KD - 3) + W 3 - f ) - 2«f - 2) - m + ^ - ^ - i)) . 

(387) 

Note, that the logarithms in / appear as a result of integration over longitudinal momenta. 
For the two-loop including the quark loop we obtain the following result [59]: 

- N -^f r (2 - / ( 2(^)- 4 - AH • 



(2tt)^- 2 v 2 ' T(D) J gf (t - g?) 2 V V V 



6 Conclusion 

As it was shown above, the theory of the high-energy scattering in QCD is based on 
the BFKL equation summing the leading logarithmic terms in the perturbation theory. 
Next-to-leading terms were calculated using the /c^-factorization only for the anomalous 
dimensions describing the transition of quarks to gluons [60]. Above we expressed the 
next-to-leading real contribution to the BFKL kernel in terms of squares of the amplitudes 
for the production of gluons and quarks with definite helicities. All infrared divergencies 
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were extracted from these expressions in an explicit form and are regularized in the IP- 
dimensional space. These divergencies cancel with the analogous divergencies from the 
virtual corrections to the BFKL equation which were also written above. It is known [55], 
that the total next-to-leading corrections to the integral kernel can be expressed in terms 
of the dilogarithm integrals and will be available soon. The effective field theory for the 
quasi-multi-Regge processes reviewed above can be used in particular for the unitarization 
program. Because the effective action is expressed in terms of the Wilson contour integrals, 
one can attempt to include into consideration also some non-perturbative effects. The 
effective field theory can be generalized to the case of the high energy behaviour of the 
amplitudes with non-singlet t-channel exchanges including one with the baryon quantum 
numbers and to the small-rr asymptotics of the structure functions for the polarized ep 
deep-inelastic scattering [61]. 
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Appendix 



There are the following helpful relations for hypergeometric functions: 
F(2h) 



F(h, h,2h;x) 



(_ i n (i _ x ) + 2^(1) - 2^{d) -d a -d b - 2d c ) F(h, h,l;l-x), 



x h {d a + d b + 2d c ) F{h, h,l]l-x) = x l ~ h {d a + d b + 2d c ) F{1 -h,l-h,l;l-x), 

x h F{h, h,l;l-x) = x x ~ h F(l - h, 1 - h, 1; 1 - x) , 

where the derivatives d a , d b and d c act on the corresponding arguments of the function 
F (122) with the subsequent substitution of a, b and c by their given values. Using these 
relations one can continue analytically G un {p~i, p|; frf,p^>) to the region near x — 1: 



T(h)T(h) 



x h x* h G un {pl,P2;p?,p2 ! ) 



T{2h)T{2h) 



= 2tt cot(7i7i) (In 1 1 - x\ 2 - 4^(1) + d a + d b + 2d c )F(h, h,l;l- x) F{h, h,l,l-x) 

+4 (i){h)i)(h) - V>(1 - h)ijj{l - h)) F{h, h, 1; x)F{h, h, l;x*). 
Analogously according to the relations 
T{2h) 



F{h, h,2h;x) = 



-x) h (ln(-x) + 2^(1) - 2^{d) -da-db- 2d c ) F{h, l-h,l;l), 
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(d a + d b + 2d c ) F{h, 1 - h, 1; 1/x) = (9 a + <9 fe + 2<9 C ) F(l - fc, fc, 1; 1/x) , 
F(/i, 1 - /i, 1; 1/x) = F(l - /i, /i, 1; 1/x) 
one can continue G un (pi, pZ; frf, p$) at large x: 

2 

G ra (p?,p2;p^,pJ) 



r(/i)r(/i) 



V{2h)V{2h) 

-2tt cot(7r/i) (in |x| 2 + 4^(1) -d a -d b - 2d c ) F(h, l-h,l; 1/x) F(h, 1 - M; 1/x) 
+4 U(h)ip(h) - - h)ip(l - h)) F(h, l-h,l; 1/x) F(h, l-h,l; 1/x) . 
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